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Abstract. We consider globally hyperbolic flat spacetimes in 2+1 and 3+1 dimensions, where a uniform light 
signal is emitted on the r-level surface of the cosmological time for r — ^ 0. We show that the intensity of this 
signal, as perceived by a fixed observer, is a well-defined, bounded function which is generally not continuous. 
This defines a purely classical model with anisotropic background radiation that contains information about 
initial singularity of the spacetime. In dimension 2 + 1, we show that this observed intensity function is 
stable under suitable perturbations of the spacetime, and that, under certain conditions, it contains sufficient 
information to recover its geometry and topology. We compute an approximation of this intensity function in a 
few simple examples. 
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1. Introduction 

1.1. Motivation. There is considerable interest in the cosmic background radiation as an indicator of the 
history and structure of the universe. In particular, the anisotropy of the cosmic background radiation contains 
valuable information. It is explained by quantum fluctuations at the early time in the history of the universe, 
whose classical remnants became visible when the universe became transparent to electromagnetic radiation after 
decoupling. Background radiation has also been used to determine the topology of the universe [TM [51 IT51 ITS] . 
This requires a theoretical model of a spacetime, in which light signals are emitted near the initial singularity 
and received by an observer. 

Independently of the physical origin of the background radiation, the question to what degree the geometry 
of a spacetime can be reconstructed from a light signal emitted near its initial singularity is of interest both 
from the mathematics and the physics perspective. However, little is known about the behaviour and properties 
of such light signals, even for simple examples such as constant curvature spacetimes in 3+1 dimensions or 
lower-dimensional models. 

In this article, we investigate the properties of such a background radiation for a class of simple examples, 
namely flat globally hyperbolic spacetimes in 3+1 and 2+1 dimensions, in which a uniform light signal is emitted 
from a hypersurface of constant cosmological time e where e — >■ 0. The observer who receives the signal at a 
cosmological time T > e is modeled by a point p G A/ of cosmological time T and a unit, future-oriented timelike 
vector w, which specifies his velocity. This leads to an intensity function on, respectively, and S^, which 
depends on both the spacetime M and the observer (p, v). 

Although this is a purely classical model without quantum fluctuations, it turns out that the resulting 
intensity functions have rich and subtle properties and contain interesting information about the underlying 
spacetimes. In particular, the intensity function exhibits anisotropics, which are due entirely to the classical 
geometry of the spacetime M or, more specifically, the initial singularity of the universal cover of M, which is 
a domain of dependence in Minkowski space. 

If the spacetime M is conformally static, i. e. characterized by a linear holonomy representation, then the 
associated intensity function is isotropic, and contains no relevant information on M. However, as soon as the 
holonomy representation of M has a non-trivial translation component, which corresponds to a universe whose 
geometry changes with the cosmological time, the intensity function contains essential information about the 
underlying spacetime. Under certain conditions, this information allows the observer to recover the geometry 
and topology of M as well as his motion relative to the initial singularity. 

The situation is particularly accessible in 2+1 dimensions, where any globally hyperbolic vacuum solution 
Einstein's equations is a fiat globally hyperbolic spacetime of the type considered in the paper and where the 
classification results by Mess allow one to explicitly construct examples of such spacetimes. Although 2+1- 
dimensional gravity is not directly related to cosmological observations in 3+1 dimensions, it plays an essential 
role as a model in quantum gravity (see [S] and references therein). 

In particular, the results in this paper could have important applications in relating the diffeomorphism 
invariant observables of 2+1-gravity to measurements with a direct physical interpretation. The former are 
given in terms of the holonomy representations which characterize the initial singularity of the spacetime and 
serve as the fundamental building blocks in the quantization of the theory. Extracting their values from the 
measurements of an intensity function thus provides a geometrical interpretation and could be used to model 
measurements in a quantum theory of 2+1-gravity. The link between the observables of 2+1 gravity and 
measurements of light signals has been explored to some degree in |20] . but measurements of background 
radiation provides a more realistic model. 

1.2. The intensity function. After recalling the relevant background material on flat maximally globally 
hyperbolic spacetimes and their description in terms of domains of dependence in Section [51 we then introduce 
the rescaled intensity function of a domain of dependence M in Section [3l The rescaled intensity function 
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is defined in Section 13.11 It is given in terms of tlie limit e — >■ of a uniform light signal emitted from the 
hypersurface of cosmological time T — e of M and received by a free-falling observer in the spacetime. Section 
13.21 and Section [3.31 contain an explicit description of the intensity function for the two basic examples, namely 
the future of a point and the future of a spacelike line. 

These two basic examples are the essential building blocks in the analysis of the intensity function for a 
general domain M, whose properties are investigated in Section 13.41 The essential result is Proposition 13.11 
which asserts that the intensity function is well-defined and locally bounded. It should be stressed that even at 
this point, the mathematical analysis of the intensity function is not as simple as it may appear at first sight, 
and some care is needed. The situation is simpler for domains M whose initial singularity is closed, which are 
investigated in Section 13.51 In this case, the associated intensity function is continuous. Note, however, that 
many relevant examples are not of this type. 

Section r3 .61 analyses the properties of the intensity functions for generic domains in 2-1-1 dimensions. This case 
is more accessible than its 3-|-l-dimensional counterpart due to a simple description of domains of dependence, 
discovered by Mess |19) . in terms of a measured lamination on the hyperbolic plane. Domains of dependence, 
which are the universal covers of globally hyperbolic flat spacetimes, are obtained from measured geodesic 
laminations on closed hyperbolic surfaces. By using these results we prove (Proposition 13. 6|) that the intensity 
function is lower semi-continuous, and that its discontinuity set is meagre. 

1.3. Stability. In Section |4] we investigate the stability properties of the intensity function. We analyze the 
variation of the intensity function under small deformations of the domain of dependence M and of the ob- 
server. Stability of the intensity function at least with respect to small changes of the observer are a minimum 
requirement for assigning any physical meaning to it. 

Again, this question turns out to be more subtle than it appears and some care is needed in the analysis. 
This is illustrated in Section 14.11 where we show in a very simple example that if a domain of dependence M 
is the limit (in the HausdorfF sense) of a decreasing sequence of finite domains Af„ (domains which are the 
intersection of the futures of a finite set of lightlike planes) the intensity of M does not necessarily coincide with 
the limit of the intensities of the domains Af„ . 

With this example in mind, we introduce in Section r4.2l a notion of domain of dependence with a flat boundary. 
Important examples of this arc finite domains, which always have flat boundary, and universal covers of globally 
hyperbolic flat spacetimes in 2-1-1 dimensions fProposition 14.22")) . We prove (Theorem 14. 131) that if a sequence 
of domains with flat boundary converges to a domain with flat boundary, then the limit of the intensities is the 
intensity of the limit. 

If a domain D does not have fiat boundary, and if (Z?„)„gN is sequence of domains with flat boundary 
converging to _D, then the sequences of the intensities t„ of I?„ always converges to a limit intensity Lum- This 
limit intensity is not necessarily equal to the intensity t of D, but it is independent of the sequence (-D„)„gN- 
We prove fTheorem l4.14[) that intensity t is at least equal to the limit intensity ium, and at most equal to dLum, 
where d is the dimension of the spacetime. In particular, this shows that the example of Section 14.11 exhibits 
the worst possible behavior with respect to this limit, as the ratio of the intensity of the domain to its limit 
intensity is the largest possible in dimension 2 + 1. 

1.4. Recovering the spacetime geometry and topology. In Section [S] we turn to the question of re- 
constructing the geometry and topology of a globally hyperbolic flat spacetime M from the intensity of the 
background radiation as seen by an observer. We investigate this question in 2+1 dimensions, and there are 
two basic remarks regarding the general situation. 

• Reconstructing the geometry or topology of the spacetime from the observed intensity function is only 
possible for observers with a sufficiently large cosmological time. If the observer is too close to the initial 
singularity, then he might see only a small part of M and could then infer little from the background 
radiation she observes. 

• The observer can only determine parts of the initial singularity of the universal cover M of M from the 
measured intensity function. Therefore, there is no way for the observer to be sure, at any given time, 
that what he observes is really the topology of M. It could happen that M is "almost" a finite cover of 
a globally hyperbolic flat spacetime M' , with only a tiny difference in a part not "seen" by the observer. 
In this case the observer could only conclude that the spacetime is either M' or one of its finite covers. 

To obtain results we make a (presumably) technical assumption, which simplifies the situation to some extend. 
We only consider spacetimes obtained by "grafting" a hyperbolic surface along a rational measured lamination, 
that is, a measured lamination with support on a finite set of simple closed curves. Under those hypothesis, 
we prove fProposition 15.81) that the observer can reconstruct the part of the lamination corresponding to the 
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part of the initial singularity that intersects his past. We also prove (Proposition 15.13")) that the observer can 
reconstruct the whole geometry and topology of the spacetime in finite eigentime up to the above-mentioned 
problem with finite covers. 

1.5. Computations for examples. In Section [HI we present explicit computations for the intensity function 
seen by an observer for three different globally hyperbolic flat 2-|-l-dimensional spacetimes. The spacetimes 
considered are chosen for their simplicity. Two are obtained by grafting a hyperbolic surface along a rational 
measured lamination, the third by grafting along an irrational lamination. For each of those spacetimes, we 
provide pictures of the intensity function as seen by an observer located at different points in the spacetime. 
This allows one to observe explicitly the variation of the intensity function depending on the cosmological time. 

In dimension 3 + 1, we only consider one example, described in Section [71 This is due to the fact that 
3+1-dimensional globally hyperbolic flat spacetimes are much more difficult to construct than their 2 + 1- 
dimensional counterparts. In both cases they are associated to first-order deformations of the flat conformal 
structure underlying a hyperbolic manifold. However, hyperbolic manifolds are flexible in dimension 2, while 
they are rigid in dimension 3,. Consequently, it becomes more difficult to find an adequate deformation cocycle 
in dimension 3 + 1. The example considered in Section [7. II is due to Apanasov [2, , and it has the relatively rare 
property of admitting several distinct deformation cocycles. We provide some pictures of the intensity seen by 
an observer in a spacetime constructed from this example. 

1.6. Possible extensions. In this article, we consider only flat globally hyperbolic spacetimes. However, it 
should be possible to perform a similar analysis for globally hyperbolic de Sitter or anti-de Sitter spacetimes, 
which have a similar structure, at least with respect to the geometry of their initial singularity. 



2. Globally hyperbolic Minkowski space-times 

2.1. Minkowski space and domain of dependences. In this section, we collect some properties of 
Minkowski space and refer the reader to [ini [Z] for details. Minkowski space in n + 1 dimensions, denoted 
R^'" in the following, is the manifold M."^^ equipped with the flat Lorentzian form 77 = —dxQ + dxi + . . . + dx^, 
often referred to as Minkowski metric. 

Isometry group. Isometrics of Minkowski space are afflne transformations of R""'""'^ whose linear part preserves 
the Minkowski metric. We denote by 0(1, n) the group of linear transformations of E""*"^ which preserve the 
Minkowski metric (Lorentz group in n + 1 dimensions) and by Isom(n, 1) the group of isometries of Minkowski 
space (Poincare group in n + 1 dimensions). 0{n, 1) is a n(n + l)/2-dimensional Lie group with 4 connected 
components, and we denote by S0^(jl,1) its identity component, which contains linear orthochronous trans- 
formations with positive determinant. The dimension of Isom(n, 1) is n{n + l)/2 + (n + 1) = (n + l)(n + 2)/2 
and for n > 3 has four connected components. The identity component, denoted Isomo(l,n), contains the 
transformations that preserve both the orientation and the time orientation. 

Flat spacetimes. It is well-known that every flat spacetime is locally modeled on Minkowski space. For globally 
hyperbolic flat spacetimes a more precise result holds (see [191 II])- For every flat spacetime M with a closed 
Cauchy surface, there is a discrete group of isometries F C Isomo(n, 1) and a convex domain D C M^'" such that 
D is F- invariant and M embeds into the quotient D/T. The domain is a domain of dependence, in the sense 
that it is the intersection of the futures of one or more lightlike planes. Domains of dependence play an essential 
role in this paper, and will be described in more detail below. The quotient space D/T is called a maximal 
globally hyperbolic flat space-time with compact Cauchy surface, for which we use the acronym MGHCF. 

Hyperbolic representations. The unit timelike vectors in M.^''^ form a smooth hypersurface, H C M.^'", which 
contains two connected components: the component that contains future oriented unit vectors, and H~ that 
contains past oriented unit vectors. Both iJ+ and H~ are achronal spacelike smooth surfaces. The Minkowski 
metric induces a Riemannian metric of constant curvature —1 on H'^ and , and equipped with this metric 
and H~ are isometric to the n-dimensional hyperbolic space M". The group SO^{l,n) acts by isometries 
on , and it is identifled with the identity component of the isometry group of . Every geodesic of _ff + is 
given as the intersection of with a timelike linear 2-dimensional plane. 
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2.2. Domains of dependence. A domain of dependence (called regular domain in [7|) D is a convex domain 
of M^'" that is given as the intersection of the future (or the past) of a number of lightlike n-planes. We will 
exclude two limit cases: the whole space and the future of a single lightlike n-plane. In other words, we require 
that R^'" \ D contains at least two non parallel lightlike n-planes. 

Simple examples of domains of dependence are the future of a point, or the future of a spacelike (n— l)-plane, 
whereas the future of a spacelike n-plane is not a domain of dependence. Geometrically interesting examples 
are the universal covers of maximal globally hyperbolic flat manifolds with compact Cauchy surface (GHMFC 
manifolds). Figure [T] shows two of those more complex examples, corresponding to the domains of dependence 
described in sections 16.1.21 and 16.1.31 




Figure 1 . Two examples of domains of dependence 

Let us recall that for any Lorentzian manifold Af , the cosmological time is a function t : M — > (0, +oo] whose 
value at a point p G Af is the suprcmum of the length of causal curves in M ending at p: 

t{p) — sup{i?(c)|c is a casual curve ending at p} . 

One of the main features of domains of dependence is that their cosmological time is a regular function: this 
means that t is finite-valued and C^'^. In fact, if Z? is a domain of dependence and p G D, there is a unique 
point r = r{p) G dD n I~(p) such that t{p) — \p — r\. Level surfaces Ha = r~^(a) of the cosmological time are 
spacelike Cauchy surfaces, and their normal vector at a point p E Ha is the vector p — r{p). 

Example 2.1. 

• If is the future of a point ro, then r{p) = tq for all p G D, and the cosmological time t{p) coincides 
with the distance of p from rg. In this case, the cosmological time is smooth (real analytic in fact), and 
the induced metric on the level surface Ha has constant curvature —l/a^. 

• If Z? is the future of a spacelike affine plane Iq of dimension k < n—1 then D is a domain of dependence. 
For p G D, r{p) is the intersection point of Iq with the affine subspace orthogonal to Iq passing through 
p. Also in this case r is smooth. The level surface Ha are isometric to M''' x H"^''' (if n = 2 and fc = 1, 
this implies that the metric is flat). 

• If £) C M"'^ is the future of a spacelike segment [po,Pi], then D is divided into three regions by two 
timelike hyperplanes Po,Pi orthogonal to [poiPi] and passing, respectively, through po and pi. The 
first region is the half-space Dq bounded by Pq which does not contain pi , the second is the half-space 
Di bounded by Pi which does not contain po, and the third is the intersection V of the other two 
half-spaces bounded by Pq and by Pi . 

For p £ Do, one has r{p) ~ po, for p e Di r{p) = pi and for p £ V, r{p) is the intersection point of 
[po,Pi] with the plane orthogonal to [poiPi] that passes through p. In this case r is smooth outside the 
boundaries of the regions Dq, Di and V and is only C^'^ on their boundaries. Level surfaces are divided 
in three regions: the regions iZo(a) = Ha n Dq and Hi{a) = Ha H Di are isometric to half-spaces of 
constant curvature — 1/a^, while Ba = Ha n F is isometric to the product of the hyperbolic space of 
dimension n — 1 with an interval of length equal to \pi — po\. (For n = 2, this is a flat strip of width 
\Pi - Pol)- 

These examples are illustrated in Figure [51 
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2.3. The boundary of a domain of dependence and the initial singularity. In this section we recall 
important facts about the geometry of the boundary dD of a domain of dependence D C E^'". 

We start by summarizing a useful description of the boundary dD. If P is a spacelike [n — l)-dimensional 
plane in R^'", the orthogonal projection tt : dD P is 1-to-l, so dD can be regarded as the graph of a convex 
function u on P. Since dD is achronal, u turns out to be 1-Lipschitz. More precisely, one finds that the graph 
of a convex function on P is the boundary of a domain of dependence if and only if ||grad7i|| = 1 at each point 
where u is differentiable. 

For any point r € dD, there is at least one future directed lightlike half-line I contained in dD that passes 
through r. It is important to note that there are always points on dD from which at least two lightlike half-lines 
contained in D originate. Indeed this occurs exactly when p is in the image of the map r : Z) — > dD introduced 
in the previous section. This subset is called the initial singularitij^ of the domain D, and will be denoted by 
T. It is the smallest subset of U such that D = /+(T). 

If one regards D as the graph of a convex function u, lightlike lines in dD correspond to integral lines of the 
gradient of u, whereas the initial singularity corresponds to the set of points in which u is not differentiable. In 
many interesting cases, the shape of the initial singularity can be quite complicated. For instance, it was shown 
by Mess [H] that for the universal covering D of a generic (2-|-l)-dimensional MGHFC spacetime, the image of 
r is a dense subset of dD. On the boundary dD, we consider the pseudo-distance defined as follows: 

• given a Lipschitz arc k contained in dD, its velocity (defined a.e.) is not timelike. So we can define the 
length of k as 

£{k) = J ^{k{t),k{t))At. 

• Given ri,r2 G dD the space of Lipschitz arcs K,{ri, r2) joining them is not empty. So we can define 

rfo(ri,r2)-inf{f(fc)|fce/C(ri,r2)} . 

^Note that the use of the term initial singularity here differs from the one common in the physics literature. The set of points 
of D at which causal curves cannot be extended into the past is the entire boundary dD. 
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As the boundary of D contains lightlike segments (whose length is clearly 0), the pseudo-distance d is not a 
genuine distance, since there exist pairs of points with d{p, q) = 0. However, the following lemma shows that 
this can occur only if the images of p and q under the map r coincides. 

Lemma 2.2 ([7]). If D is a domain of dependence andp,q G D then r(p) ^ r(g) implies do(r(p), r(q)) ^ 0. 

In other words, this lemma states that the restriction of do to the initial singularity is a distance. Let us also 
remark that the topology we consider on T is the one induced by the distance do, which in general is different 
from the topology induced by Minkowski space. 

Example 2.3. 

• If Z) is the future of a point rp, the initial singularity contains only one point that can be identified with 
ro- 

• If Z? is the future of a affine subspace E of dimension k < (n — 1), then the initial singularity is isometric 
to E. 

• If D is the future of a segment in R^'^, then the initial singularity is the segment itself. 

• If C M^'^ is the intersection of three half-spaces bounded by lightlike planes, the initial singularity is 
the union of three spacelike rays starting from the intersection point of the planes. Note that in this 
case, the initial singularity is not a submanifold. In fact, generically the initial singularity does not have 
a manifold structure. The geometry of the initial singularities of domains of dependence in dimension 
2-1-1 will be studied in more depth in the next section. 

2.4. The 2 + 1-dimensional case: the Mess construction. Mess [H] discovered an efficient way to construct 
regular domains in R^'^. This construction is general in the sense that every regular domain can be constructed 
in this way. It also has the major advantage that the geometrical features of the initial singularities are readily 
apparent. 

We will describe the Mess construction in simple cases, namely for domains obtained by grafting along 
weighted multicurves. These simple cases are dense, in the sense that every domain of dependence can be 
approximated by domains of dependence obtained in this way. 

Let us start from a collection of disjoint geodesies of L = li U . . . U Ik, and a collection of positive 
numbers ai, ....a/j. Every geodesic is given as the intersection of with a timelike linear plane Pi, . . . , P^. 
The planes Pi disconnect iZ+ into a collection of regions Di, . . . ,Dh. (Note that each of them is the cone on 
some component of \ L). 

For each plane Pi, let Vi be the vector in M}''^ characterized by the following conditions: 

• it is orthogonal to Pi with respect to the Minkowski metric (in particular it is spacelike); 

• its norm is equal to a^; 

• it points to the component of R^'^ that does not contain Di. 

Now for any region Dj take the sum of all vectors Vi associated with to planes Pi that separate Dj from Di : 

i:Pi separates Di from Dj 

Translating each region Dj by the vector Wj yields a collection of disjoint domains D'l, . . . , D'^^ which are convex 
cones with vertices at Wj. In particular, note that if Z?ji is adjacent to Dj^, then Wj-^ ~Wj2 9- vector orthogonal 
to the plane Pi separating Dj-^ from Dj^ and is of norm ai. 

In order to connect the domains D'^ we consider the domains Vi obtained as follows. If Dj-^ and Dj,^ are 
adjacent along Pi, then Vi is the region of the future of the segment Si = [wj-^ , Wj^] bounded by the two timelike 
planes orthogonal to the segment Si through its end-points. 

It turns out that _D = IJ Z?^ U IJ is a domain of dependence. The map r can be easily defined on each piece: 
r sends all points of D'^ onto Wj , while it sends points of Vi to the segment Si = [wj^ , Wj^ ] . The level surface Ha 
can be decomposed into different regions: the regions Ha H D[, which have constant curvature — and the 
regions Ha H Vi , which are Euclidean strips of width . 

The initial singularity is then given as the union of the line segments Si and the vertices Wj. In particular, 
it is a graph with a vertex for every region of ZZ+ \ L. Two vertices Wj-^ and are connected by one edge if 
and only if the corresponding regions are adjacent. Combinatorially the singularity is a tree, that is, a graph 
which does not contain any closed loop. Notice that the length of each segment Si is precisely ai. 

Although we summarized this construction for a finite number of geodesies, it works analogously also when 
L is an infinite family of disjoint geodesies that is locally finite (i.e. every compact subset of ZZ+ meet only a 
finite number of li). 
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2.5. The equivariant construction. Using tlie construction from the previous subsection, one can construct 
the universal coverings of MGHFC spacetimes different from I"'"(0) as follows. Take a hyperbolic surface F 
and consider the metric universal covering tt : iJ+ — ?► F and covering group F < S'0+(l,2). Consider on 
F a disjoint collection of simple closed geodesies ci . . . Cfc and positive numbers ai . . .a^. Then the preimage 
L = 7r^^(ci U . . . U Cfc) is a union of infinitely many disjoint geodesies. The weight of each geodesic k C L is the 
number corresponding to n(li). As above, the geodesies k correspond to planes that cut I'*"(0) into infinitely 
many pieces Dj. By the invariance of L under the action of T, elements of F permute the regions Dj. 

The construction explained in the previous subsection then produces a domain D, and Mess showed that 
there is an affine deformation F' of F, so that D is F'-invariant and the quotient is a MGHFC spacetime. 
Namely any 7 e F is changed by adding a translation part of vector ^(7) which is the sum of all vectors Wi 
corresponding to the planes Pi disconnecting Di from jDi. 

Remark 2.4. In the example above, it can be seen that each Dj bounds infinitely many planes Pi. This implies 
that the vertex in the initial singularity corresponding to Dj is the end-point of infinitely many edges, or 
equivalently has infinite valence. 

In the examples illustrated in the previous section it turns out that the initial singularity of domains of 
dependence in R^'^ is always a graph, and in fact a tree (possibly with vertices of infinite valence). In fact, 
there are more complicated examples in which the initial singularity does not have a simple graph structure, 
but it is always a real tree according to the following definition. 

Definition 2.5. A metric space (T, d) is a real tree if for every p,q E T there is a unique arc k C T joining 
them. Moreover k is the image of an isometric immersion I ^ T where / is an interval of length equal to d{p, q). 

Real trees are generalizations of the usual trees (which, by contrast, are often called simplicial trees). The 
domains of dependence whose singularity is a simplicial tree are exactly those constructed in the previous 
section [Sl|3]. In particular, every domain of dependence with simplicial tree as initial singularity is determined 
by a simplicial measured geodesic lamination of H'^ , which, by definition, is a locally finite union L of disjoint 
geodesies li, each equipped with a weight > 0. 

Proposition 2.6. ^ If D is a domain of dependence in R^'^ then its initial singularity T is a real tree. 
Moreover, the vertices of T are those points in dD at which at least three lightlike segments in dD originate. 

Given a point r let Dr be the convex hull in Minkowski space of the lightlike lines contained in dD 

which start at r. Then Dr is a convex subset of I"*"(?'). Notice that the dimension of Dr is 3 if and only if r 
is a vertex, otherwise Dr is the intersection between D and the timelike plane containing the two lightlike rays 
starting at r. 

If Tr is the translation which send r to 0, we denote by J> the intersection of with Tr{Dr). Note that J> 
can be interpreted as the set of unit normals of the support planes of D at r. A number of consequences follow 
directly. 

• If r is a vertex then J>. is a region of bounded by disjoint geodesies. 

• If r is not a vertex then is a complete geodesic. 

• If r 7^ s then and J^, have disjoint interiors. Jv and J-g can be disjoint, they can coincide if they are 
both lines, or they can meet along a boundary component. 

In particular the set L = [j^^^^ ^^^.^^^ dTr U IJ^ is not a vertex -^r is a union of disjoint geodesies that are called 
the leaves of L. In general, the set L can be quite complicated. The intersection of a geodesic arc in i7+ and L 
can be uncountable (and sometimes a Cantor set). 

The simplest ease is when the singularity is a tree. In this case, the set L is the union of isolated geodesies: 
any compact subset of meets only a finite number of leaves. In this case, it is also evident that components 
of \ L corresponds to vertices of T, whereas each leaf of L corresponds to an edge of T. 

In addition to L we can construct a transverse measure that is the assignment of a non-negative number for 
any arc transverse to the leaves of L which verify some additivity conditions, see e.g. (Sj. If A: is an arc on 
i7+that joins two points in i7+ \ L and meets each leaf at most once (for instance if fc is a geodesic segment 
that is not contained in any leaf), we define fi{k) ~ do(^Oj ^1) where rg and ri are the points on T such that the 
end-points of k are contained in J>(, and J>j^ . 

If the lamination is locally finite, for each leaf I there is a number a{l) that coincides with the measure of any 
are k transversely meeting only I. Any transverse arc k can be subdivided into a finite number of arcs /ci , . . . , kp 
such that each ki meets every leaf at most once. So we can define ^(fc) = J^l^i^i)- 

Mess [in] showed that the data (L, /j) determines D up to translation. In the simple ease where the lamination 
is locally finite, the construction of D from (L, p) is the one summarized in Section [2.41 
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Remark 2.7. In dimension n + 1 > 4, it is no longer true that the initial singularity is a tree. In fact the 
geometry of the initial singularity is still not understood. In [7] a description of the singularity is given in some 
special cases. 

2.6. Holonomies of domains of dependence and hyperbolic structures. In Section [^?T1 we summarized 
the construction which assigns a domain of dependence to each flat GHMCF manifold. There is also a deep 
relation between holonomies of flat Lorentzian manifolds and first-order deformations of holonomy representa- 
tions of hyperbolic manifolds, which was already used in dimension 2-1-1 for instance in |12| . In the following, 
we summarize this relation, which behaves somewhat differently in dimension 2-1-1 and in higher dimension. 

Dimension 2+1. In this subsection, we recall how a flat (2-|-l)-dimensional GHMCF manifold can be obtained 
from a point in Teichmiiller space together with a deformation 1-cocycle. For this, note that K^'^ can be identified 
with the Lie algebra sl{2, M) with its Killing metric. The canonical action of SO{2, 1) on M."^'^ corresponds to the 
adjoint action of SL{2, R) on sl{2, M). For each representation of 7ri5 in SO{2, 1), it determines a vector bundle 
over S with fiber M^'^, which corresponds to the sZ(2, R)-bundle over S defined by the adjoint representation. 

Proposition 2.8. |19j Let M be a flat GHMCF manifold homeomorphic to S xR, where S is a closed surface of 
genus at least 2, and let h : ttiS — > Isom(2, 1) be its holonomy representation. Then h decomposes in Isom(2, 1) = 
50(2,1) K M.'^'^ as h = {p,t) where p : ttiS 5*0(2,1) has maximal Euler number, and r : ttiS — >■ sZ(2,R) 
is a 1-cocycle for p. Conversely, any couple {p,t) where p : ttiS — > 50(2,1) has maximal Euler number and 
T : TTiS — ?> s/(2,R) is a 1-cocycle for p defines a representation of ttiS in Isom(2, 1) which is the holonomy 
representation of a flat GHMCF manifold. 

One way to obtain a 1-cocycle is by considering first-order deformations of a surface group representation 
in 50(2,1). This is summarized in the following proposition, which allows one to construct the holonomy 
representation of a fiat GHMCF is as a first-order deformation of the holonomy representation of a hyperbolic 
metric on a surface. 

Proposition 2.9. Let {pt)te[o,i] be a smooth one-parameter family of morphisms from 7ri(5) to P5L(2,R). 
Then the map r — pQ^{dp/dt)t=o from 7ri(5) to s/(2,R) is a 1-cocycle for pQ. 

GHMCF spacetimes as first-order deformations in higher dimension. We will now consider the construction of 
flat GHMCF spacetimes as first-order deformations in dimension n -I- 1 > 3. For this, we consider a closed, 
orientable, hyperbolic n-dimcnsional manifold M , with fundamental group F. The holonomy representation of 
M is a homomorphism po : T ^ SO{n, 1). It is rigid by Mostow's theorem, and also infinitesimally rigid, in the 
sense that any deformation cocycle for po vanishes. 

However, M can be considered as a totally geodesic hypersurface in a complete, non-compact hyperbolic 
manifold N of dimension n -\- 1. This corresponds to extending pQ to a representation p : F — > 50o(»t. + 1,1) 
with image in SOo{n, 1) C SOo{n + 1, 1). Now consider a deformation {pt)te[Q.i] of P- in dimension 2 + 1 
one obtains the map 

pi :=p(0)-V'(0) :r^o(n + l,l) , 
which is a deformation cocycle for p. Moreover, there is an orthogonal decomposition o(n + l, 1) = o(n, 1)®R"'"'^, 
and we can decompose pi along this direct sum. The component in o(n, 1) is a deformation cocycle for po, so 
it vanishes by the infinitesimal rigidity of po, and thus pi determines a R"' ^-valued cocycle. 

This cocycle then determines a GHMCF spacetime, with holonomy representation (po,pi) considered as 
a homomorphism from F to Isom(R"'^) — SOo{n,l) k R"'^. Moreover, the holonomy representations of all 
GHMCF spacetimes can be obtained in this way, see [151 [7] . 

Proposition 2.10. The GHMCF spacetimes for which the linear part of the holonomy is equal to po are in 
one-to-one correspondence with the deformations of p. 

There is another geometrical interpretation of the deformation cocycle p, namely as a first-order deformation 
of the flat conformal structure on M underlying its hyperbolic metric. Indeed, in the situation described above, 
where M is considered as a totally geodesic submanifold of N, the conformal structure at infinity of N remains 
conformally fiat under a deformation. Conversely, any first-order deformation of the conformally flat structure 
on M determines a deformation of its developing map in 5" and, by taking the convex hull of the complement 
of its image, one obtains a first-order deformation of N. 

One direct consequence Proposition 12 . 101 is that it is much more difficult to construct examples of GHMCF 
spacetimes in higher dimensions than in dimension 2 + 1. An (n+l)-dimensional GHMCF spacetime is uniquely 
determined by a closed n-dimensional hyperbolic manifold along with a R"'^-valued deformation cocycle. For 
n = 2, the latter corresponds to a tangent vector to the Teichmiiller space for a surface 5 of given genus g, 
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which imphcs that the GHMCF spacetimes homeomorphic to 5 x R form a manifold of dimension 12g— 12. For 
n > 3, finding a GHMCF spacetime homeomorphic to M x R is more difficult. Any closed manifold M admits 
at most one hyperbolic metric g by Mostow's rigidity theorem. Finding a deformation cocycle is equivalent to 
finding a first-order deformation of the warped product hyperbolic metric dt^ + cosh^ (t)g on M x M. For many 
choices of {M,g), such a deformation cocycle does not exist. However there are also many examples where 
{M,g) does admit a M"'^-valued cocycle. 

• This occurs whenever {M, g) contains a closed, totally geodesic submanifold, and the cocycle corresponds 
to "bending" along this totally geodesic manifold, see [T71 [T^ . There are many (arithmetic) examples 
of closed hyperbolic manifolds (in any dimension) containing a closed totally geodesic surface. 

• Other examples of deformation cocycles can be found in specific cases, see e.g. [TBI [^. 

In Section [7] we investigate the examples constructed by Apanasov in [2] and show how the intensity function 
encodes information on the holonomy representation and hence on the topology of the spacetime. 

2.7. Reconstructing a domain of dependence from its holonomy representation. In Section [S] and [7] 
we compute domains of dependence which are universal covers of MGHCF spacetimes. This requires a practical 
way of reconstructing (to a good approximation) the shape of a domain from the holonomy representation of 
the MGHCF spacetimes. For this, we use another characterization of those domains due to Barbot ji. 

We consider a MGHCF spacetime M of dimension n + 1 with fundamental group F. As explained in Section 
12. 1[ the universal cover of M can be identified isometrically with a future-complete domain of dependence 
D C K^'". The fundamental group F acts isometrically on D with a quotient D/F isometric to M. Moreover 
all elements of F except the unit element act on R^'" as loxodromic elements. 

Definition 2.11. For g E T, we denote by Dg the set of points x E M^^" such that g^{x) — x is spacelike for all 

This is a simpler version of the definition at the beginning of Section 7 in [3] , but both definitions are equivalent 
in our case because the linear part of each nontrivial element (/ G F is loxodromic. In (2-f 1) dimensions, it 
is easy to give a more explicit description of the set Dg. If the linear part of g is loxodromic, then there is a 
unique spacelike line Ig in M^'^ which is invariant under the action of g. It is proved in [5^ that the set Dg is 
then equal to the union of the past and the future of Ig. 

Proposition 2.12 (Barbot [3 ). The domain of dependence D is one of the two connected components of 
^gerDg. 

A proof can be found — in a more general setting — in Barbot's work [3], see Section 7 for the definitions 
and Section 10 for the statements corresponding to Proposition l2.12l 

3. Light emitted by the initial singularity 

3.1. Definitions. We consider a domain of dependence M in (n-f l)-dimensional Minkowski space, as described 
in the previous section. 

An observer in free fall in M is characterized by his worldline, which is a future-oriented timelike geodesic in 
M. This geodesic is specified by the choice of a point p £ M and a future-directed timelike unit vector v £ H", 
where we use the identification of H" with the set of future directed timelike unit vectors iJ+ from section [5J 
The point p £ M corresponds to a given event on the worldline of the observer, while the vector v is his velocity 
unit vector. 

We consider a uniform light signal emitted near the initial singularity of M which is received by the observer 
at the point p £ M . The quantity measured by the observer is the intensity of this light signal, which depends 
on the observer and on the direction in which the light is observed. We can construct this quantity as follows. 
The space of lightlike rays arriving at p can be identified with the set of unit spacelike vectors orthogonal to 
V and hence with T^M", which corresponds to the (n-l)-dimensional sphere 5""^. We associate to each unit 
vector u £ T^M" the ray through p with the direction given by the lightlike vector u — v. The basic idea is to 
define the (rescaled) intensity as a function 

which is given as the renormalized limit of the functions that measure the intensity of the light emitted from 
the surface surface iJ^ of constant cosmological time e: 

Pp,v{u) = limep^ (m). 
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The functions are defined by the rule 

where ^{u) is the normal of the surface at the intersection point of with the light ray p + M.{u — v): 
{q,{u)}^H,C^{p + ^{-v + u)). 

We consider first the (2 + l)-dimensional case. In this situation, the intensity function of a domain that is 
the future of a finite spacelike tree can be understood by considering two main examples. The first is a domain 
that is the future of a point, i. e. a light cone, and the second is a domain which is the future of a spacelike line. 
We first consider these two examples and then use them as the building blocks to analyze the general situation. 

3.2. Example 1: future of a point. We consider the domain of dependence D which is the future of G K^'^ 
together with an observer in D which given by a point p £ D and a future directed timelike unit vector v £ 
as shown in Figure [31 The cosmological time t of the event p is then given by 

(1) (p,p)=-r2. 
We also consider the quantity 6 defined by 

(2) {p,v) — —TCOshS, 

which is the hyperbolic distance between v and the point p/r (see Figure [21) • and measures the discrepancy 
of the observer's eigentime and the cosmological time. The cosmological time coincides with the observer's 
eigentime up to a time shift iS S = 0. We can also interpret it as the rapidity of the boost from the worldline 
of the observer to the the geodesic through p and the origin. 

For a given unit vector u e T^lP we also introduce a parameter defined by 

(3) (p, u) = rsinh , <f) e [-5, 6] . 

Geometrically, is the hyperbolic distance of the point p/r from the geodesic orthogonal to u. It becomes 
maximal when u € Tj^H^ points in the direction of p/r e and minimal when u points away from it. 

We denote by e M"*" the parameter that characterizes the intersection point qe{u) of the light ray p+M(u— w) 
with the surface ffg of constant cosmological time e. As the latter is the set of points 

H^^{xe M^'^l {x,x) = -e^}, 

the parameter is characterized uniquely as the positive solution of the equation 

(p + t,{-V + U),p + t,{-V + U)) = -€^. 

Inserting the parameters 5 and defined in ^ and ^ and solving the equation for t^, we obtain 

^2 ^2 



2T(cosh(5 + sinh(/)(u)) 

The unit normal vector ^ in the intersection point of p + M(u — v) and is given by 

= ^(P + <e(-W + M)), 

and the function „ by 

P + ie( — V + u) " 

e 

A direct computation then shows that the rescaled intensity function then takes the form 



p.v 



Pp,v[u) 



2(cosh(5 + sinh0(u)) 

Using the fact that the function (p takes values 0(w) G [—(5,(5], one finds that the maximum and minimum 
intensity are given by 

(A\ n""^^: _ ^ 5 rain _ I„-<5 

V^J rp,v ~ 2 ' "Pi" ~ 2 

These equations show that the intensity function pp^y allows one to re-construct the cosmological time t(p) of 
the observer at the reception of the light signal and the discrepancy between his eigentime and the cosmological 
time, which is given by S. Moreover, by determining the direction of the maximum, the observer can deduce 
the direction of p/r e H^. 
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Figure 3. Description of an observer for Example 1. 



3.3. Example 2: future of a spacelike line. We consider the domain of dependence D which is the future 
of a spaceUke line Re in M^'^ together with an observer given by a point p € D and a timelike future directed 
timelike unit vector v e H'^. Using the symmetry of the system under translations in the direction of the line 
and Lorentz transformations with the line as their axis, we can parameterize the data for the observer as follows 

(5) e= (0,0,1) p = T(cosh(5, sinh5, 0) w = (coshf , 0, sinh^). 

If we denote by xq € the timelike unit vector corresponding to the shortest line in R^'^ from Re to p, then ^ 
is the hyperbolic distance of ?; S from from the geodesic through xq that is orthogonal to Re. The parameter 
5 is the hyperbolic distance from xq to the projection of v onto this geodesic, as shown in Figure H) 
The rescaled intensity function is defined as in Example 1. For u S T^IS? , we have 

(6) „ (m) = lim ep^ „ {u) „ (m) = (u, v^,^., {u)) , 

where Vp^v{u) is the unit normal vector to the constant cosmological time surface at the intersection point 
= n (p + R(m — v)). As the constant cosmological time surface is of the form 

(7) H,^{xe R^'^ I (x, x) ~ {x, ef ^ -e^} , 
the intersection point q^iu) is given by the equations 

(8) qe{u) ^ p + t^{u~ v), {p + t^{u-v),p + t^{u-v)) - {p + te{u-v),e)'^ ^ -e'^. 

Using the parameterization ([Sj, in particular the identities {u — v,u ~ v) — 0, (p, e) = 0, we obtain a quadratic 
equation in 

(9) t^^{e,u-vf - 2U{p,u-v) -e^ ^0 
with solutions 

(p, u — v) zt \ (v, u — f )^ — (t^ — e2)(e, u — v)"^ 

(10) tf = - ^ 



In the limit e — > this reduces to 



(p. 7i — u) ± 'J {p,u — v)^ — T-^{e, u — v)^) 

(11) t± ^ — -, 72 ' 

(e, u ~ V) 

and the rescaled intensity function is given by 

(12) p ^ ^{y^p) - t,{l - {u - v,e){v,e)). 

To obtain a concrete parametrization for p, we parameterize the unit vector u G T^{M^) as 

(13) M = cos6l(sinhC,0,cosh^)+sin6l(0,l,0) . 
This implies 

(w,p) = — Tcosh(5cosh^ , (w,e)=sinh^, 
(p, u — v) = — TCOsh(5(cos flsinh^ — cosh^) + rsinh S sin 9 , 
(e, u~ v) — cos 6'cosh ^ — sinh ^ , 
and the expression for t± becomes 

cosh ^ — sinh ^ cos 6* =F sin 
(cosh ^ cos 6* — sinh^)2 
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If we introduce an "angle variable" defined by 

Q 

(15) tan Y = e« with 9^ € [0, tt/2] 
then we obtain 

(16) "^"^'^^"^^ 



2 sin' 



2 f e±e, 



\ 2 



Note that t±> Q and that t± diverges for — =F%- The two cases for t±, p± correspond to the intersection 
points of the ray p + E(m — v) with the two lightlike planes Q± containing Me. The relevant intersection point 
is the one that is closer to p, i.e. the one with t — niin(t±). From it follows that this is the one associated 



with tj^ if 



and the one for t- otherwise. For the associated intensity functions, we obtain 



T , J ±s zrA /coshC-sinh^cos6'Tsin6'\^\ t /sm(^\ 
p^iO)^-cosH[e^'-e^^[ cosh ^ cos ^ - sinh ^ ) = 2^°^'^ 



2\ 



sin [-2^ ) 

Clearly, p+{0) > if and only if p7)) is satisfied, and p-{0) > otherwise. The intensity is therefore given by 



V 



2 

25 



(18) pp,,{e)=max{p+{9),p^{9)) = { 




p+{9) for I ^ ) < e 

P-{9) for 



The intensity function has local maxima in 

(19) 05^"^ = ±6*5, 
where it takes the values 

(20) = ^coshee^'', 
and it vanishes if and only if 

(21) 



Sm [-2^ ) 



This corresponds to t+ = t_ or, equivalently, 

(22) {p~ Te,u- v) = 0. 

This condition is satisfied if and only if light ray p + K(u — v) intersects the line R • e, which, for each observer, 
happens for exactly two values of 9. 

The observer can therefore extract all relevant information from the function pp^y{9). He can determine the 
cosmological time r at the reception of the light signal, his position relative to the line, which is given by S, and 
his velocity relative to the line, which is given by ^. The development of the measured intensity in terms of the 
eigentime of a moving observer is given by the dependence of his cosmological time and the parameter S on his 
eigentime. For an observer with a worldline specified by p e and v e H^, his position at an eigentime t after 
the event p is given by p' ~ p + tv. This implies that the cosmological time of p' and the associated parameter 
5 are given by 

(23) tU) = Jt'^ -2(p,v) = v/t2 +i2 +rcosh5coshC , coth S{t) = coth S + t^^^^^ . 

cosho 

The time development of r with the eigentime corresponds to an overall rescaling of the intensity function. The 
time development of 5 to changing the relation between its constant and its angle-dependent part. 
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3.4. The rescaled density for a general domain. The results in the previous subsections allow one to 
construct the rescaled intensity for a domain which can be obtained from the light cone by the Mess construction 
[19] . i. e. via grafting along a weighted multicurve. To show that the rescaled density is well-defined also for the 
case of a general geodesic lamination, one has to prove that the limit lime_>.o ep^ ,y exists for these domains and 
for a general observer (j),v). As we will see, in general the problem is more subtle than it appears and some 
care is needed to pass to the limit. 

Let D be a generic domain of dependence in R"'^. Fix an observer in D by specifying a point p E D and 
a future directed timelike unit vector v. To analyze the behavior of the intensity function, it is convenient to 
express the function epp,t,(e) as the composite of two functions : TyH^ — )> i/e and ■ ^ R, where (?e(e) is 
the intersection of the light ray p -f M • {e — v) with the level surface of the cosmological time = (e), and 

L^{q) = ~e{iy^{q),v) , 

where iy{q) denotes the unit normal vector of in q. It is clear that the maps q^^ : rjH^ converge to a 

map qo : T^lP — ?> dD as e — )■ 0. The idea is now to show that the maps te : — > M converge to a function 
c : dD — 7> R as e — >■ 0. Note, however, that functions : — >■ R are defined on different domains, so we need 
to make this statement more precise. 

Let Py be the hyperplane in R"'^ orthogonal to v. By [l9l[7], the surfaces can be realized as the graphs 
of convex functions : Py M. More precisely, the level surfaces of the cosmological time can be identified 
with the set of points ~ {q ~ x + u^{x)v\x G Py}- As e — >■ 0, the functions converge to a convex function 
uq, whose graph is the boundary dD. 

Thus, there is a natural identification between Py and i/j given hy x t-^ x + Ue{x)v. In particular, we may 
consider the functions te as functions defined on P„ . The following result then shows that the intensity function 
can be defined pointwise on the boundary of any domain of dependence. 

Proposition 3.1. For a fixed x d Py the function e bt{x) is increasing. Moreover the function 

i{x) = lim ic{x) = inf i^{x) 

e— >-0 e 

is finite-valued at each point and locally hounded. 

The proof of this proposition will be based on the following technical lemma. 
Lemma 3.2. The square of the cosmological time is convex along each timelike line. 
Proof. Take r e R"'^ and consider the cosmological time function on /^(r) which is defined by 

Trip) = \/-{p-r,p-r) 

It is clear that the restriction of along every line p + R • w is a degree two polynomial function of the affine 
parameter with leading coefficient — (w, w). In particular, the function is convex along all timelike directions. 

Given a point q (z D, let r =^ r(q) be the corresponding point on the singularity. Then I^{r) C D and on 
/+(r) we have that r > r^. Moreover, the cosmological time of q and its gradient are given by T{q) = Tr{q) and 
gradT(q) = gradTr{q) = i(g - r). 

Take a timelike vector w and consider the functions f{t) — r^(q -f- tw) and g{t) = T^{q + tw). They are 
C-'^'^-functions, which coincide with their derivatives a.t t — 0. Since f{t) > g{t) we deduce that if /" exists in 
then /"(O) > ^"(O) > 0. Thus f"{t) > for all t for which /" exists. Since /' is Lipschitz, for s < t 

f'{t)-f'{s) = lj"{x)dx>0, 
and / is convex. □ 
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Proof of Proposition \3.1i Let D be a domain with cosmological time function r. It then fohows from the resuhs 
in [7], that for ah points q G H{e) one has gradT(g) = —i/^{q). This imphes — ez^e = rgradr — igrad(r^), and 
we deduce ^ 

Le{x) = -{gradr^ix + u^{x)v),v) . 

For a given point x S P^, , we can consider the restriction of to the vertical hne a; + K • 

fx{s) = r^(x + sv), 
which is a convex function of s by Lemma 13.21 As we have 

l-eix) = ^{fx)'{u,{x)), 

the monotonicity of then fohows from the monotonicity of (fx)' ■ Q 

Proposition 13. II ahows us to define the rescaled intensity function for a domain D and an observer {p,v) as 
the map 

Pip,v) ■ rjtf R>o, p{e) = P{p,v){e) = tt,(go(e)) . 
Note, however, that it is in general not true that Pe{e) — >■ p{e) as e — )■ 0, since the convergence — > t is not 
necessarily uniform. In the next section, we will investigate the regularity of i and show that the convergence 
of pe to p holds for generic observers and generic directions. 

3.5. Domains with a closed singularity in dimension 2 + 1. To analyze the convergence properties of the 
rescaled intensity functions pe, we first consider the case where the initial singularity T is a closed subset of dD. 
In this situation, the intensity function simplifies considerably. Note, however, that this condition never holds 
for the universal covering of MGHCF spacetimes, as will be proved in the next section. Nevertheless, the results 
for this case are useful to compute the intensity function for domains of dependence that are the intersection of 
a finite number of half-spaces. 

The simplifications in the case of a closed initial singularity arise from the fact that there is an extension of 
the map r on the boundary of D, based on the following geometric idea. For each point go G dD, there is a 
lightlike ray R through qq which is contained in dD. We will suppose that the lightlike ray R is maximal with 
respect to inclusion. The ray R can always be extended to infinity in the future, but it has a past endpoint 
To € dD. The ray R is unique unless go = ''0 5 which implies that the point vq is uniquely determined by qq. 

This defines a natural map ro : dD — )• dD, and it follows directly from its definition that Tq o ro — Tq. 
The image of ro is called the extended .singularity and denoted by T. It contains all points which are the past 
endpoint of a maximal lightlike ray contained in dD. In particular, the initial singularity T is contained in 
T. Note that, unless T is closed, the map ro cannot be continuous. We will see in the next section that this 
non-continuity occurs in many interesting and relevant examples. 

Proposition 3.3. Consider the function ly : dD — > K>o defined by 

Iv{qo) = {qq - ro(go),w) • 

// the singularity is closed in dD then T — T and uniformly converges to /„. As a consequence, the function 
iy = /„ is continuous. 

Proof. For any point x S P„ we denote by qc{x) the point x -\- Ue{x)v S and by r^(x) be the projection of 
qe{x) on the initial singularity. The results in [7] imply that these points are related by the following equation 

(24) q^{x) = r^{x) + ev^{x) . 

Take now any sequence x„ S Pv that converges to x and e„ — > 0. Then the associated sequence q,i{x) 
converges to qQ{x) — x + uq{x)v, and we obtain 

LcS^n) = -{qe„{Xn) - re„{Xn),v) . 

To prove that converges uniformly to ly, it is then sufficient to check that r„ = r^^(xn) converges to ro(go) 
For this, note that the sequence r„ is contained in a compact subset of the boundary dD and hence has a 
subsequence {rn^)keN which converges to a point ro € dD. By the assumption on T, the point ro is also 
contained in T. We prove in the next paragraph that ro = r{qo{x)). The uniqueness of the limit implies that 
the whole sequence r„ converges to ro. 

That ro = ro(qo) can be established as follows. The sequence of timelike vectors q„^, — r„^, converges to 
<Zo ~ ''o- This implies that qo — ro is not spacelike. As dD is an achronal surface, it must be lightlike and the 
lightlike ray R — ro + R>o • {qo — ro) is contained in dD. Since ro is on the singularity, we obtain that the ray 
R is maximal so that ro — T^{qo)- This also shows that the image of ro is contained in T. □ 
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Remark 3.4. In the general case, we cannot conclude because the limit point tq may not be on the singularity. 
However, it is always true that the ray go + K • (go ~ ''o) is contained in dD, so it is contained in the maximal 
lightlike ray through go- In other words the point tq lies on the segment [go, ro(go)]. This shows that in general 

limsup tc„ (x„) < Ivix) 

for any sequence Cn — )■ and Xn — ^ x. In particular, this implies tyi^x) ^ Iy(^x) and hence that iy is zero on the 
initial singularity. 

3.6. Generic domains in dimension 2 + 1. Although at a first sight, the hypothesis of Proposition 13.31 could 
appear to hold generally, this is not the case. Indeed, if D is the universal covering of a MGHFC spacetime 
whose holonomy representation is not linear, then the condition cannot be satisfied. 

Remark 3.5. [TO] Let D be the universal covering of of a MGHFC spacetime. Then T is never closed in dD. 

The following proposition describes the regularity properties of of the functions i,y and Iy for in general 
domains in 2+1 dimensions. The result is that at generic points, these functions are continuous and coincide. 

Proposition 3.6. The following properties hold for the functions by and Iy: 

• The function by is upper semicontinuous. 

• The set of discontinuity points of by is meagre. 

• The function Iy is upper semicontinuous. 

Proof. The first property holds since / is the suprcmum of a family of continuous functions. Moreover, as it is 
the limit of continuous functions, by a classical result of Lcbesgue, its discontinuity points form a meagre set. 

Let us prove that Iy is upper semi-continuous. For this, take a sequence of points Xn G Py that converges 
to X. Up to passing to a subsequence we may assume that \vais\xp Iy{xn) = lim/^(.T„). If lim/„(a;„) = 0, then 
clearly Iy{x) > limsup/^(a;„). 

Let us treat the case where lim/,j(a;„) > 0. The sequence of points g„ = Xn + u{xn)v € dD converges to 
g = x + u(x)v e dD. By the assumption on lim/„(a::„) we have that ro(gn) ^ g,i for n sufficiently large. Consider 
the sequence of lightlike rays i?„ containing g„ and ro(gn). Up to passing to a subsequence, we may assume 
that it converges to a lightlike ray R through g. The sequence (ro(g„))ngN converges to the past end-point of 
R. Since R is contained in dD we deduce that ri — limro(g„) is a point on the segment [ro(g), g] and we have 

lim Iy{xn) = -(q-ri,v) < (g-ro(g),i') = Iy{q) ■ 

n^-\-oo 

□ 

4. Stability of the intensity function 

In this section we investigate the stability of the intensity function. Given a sequence of domains of depen- 
dence Dn that converges to D and an observer (p, v) in D, then (p, v) is also an observer in D„ for n sufficiently 
large. Let p" be the intensity function for Z3„ as seen by the observer (p, v) and p the associated intensity 
function for D. We will investigate under which conditions the intensity functions p" converge to the intensity 
function p. 

We will first show by a counterexample in subsection 14 . II that in general there is no convergence even in the 
weak sense. However, we identify a subclass of domains of dependence, called domains of dependence with a 
fiat boundary, which includes the interesting examples. In subsection 14.21 we prove that the convergence holds 
for these domains. In subsection 14.51 we will then prove that universal coverings of MGHFC spacetimes in 
dimension 2+1 are contained in this class. As in the previous subsection, it is advantageous to work with the 
intensity functions by and i" introduced there. 

4.1. An example. We fix coordinates xo,xi,X2 on R^'^, so that the Minkowski metric takes the form —dxQ + 
dxi + dx2 and consider the timelike vector v — (1,0,0). Let P be the horizontal plane at height equal to 1. 
Then the intersection of P with the cone I'^{0) is a circle C of radius 1. Let Ck be the regular polygon with k 
edges tangent to C. Clearly, Ck converges to C in the Hausdorff sense as k ^ +oo. 

Now observe that for each edge of Ck the plane that contains and this edge is lightlike, since it is tangent 
to I~^{0). Let Dk be the intersection of the future of the k lightlike planes containing and the edges of Ck- 
Then Dk is a domain of dependence and converges to D for fc — ^ oo on compact subsets. 

We denote by ijj and by, respectively, the intensity function of Dk and D with respect to the observer {p,v), 
and we will prove that i^' does not converge to by on a set of positive measure. Regarding L^-functions as 
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continuous functionals on the set of continuous functions with compact support, we will show that weakly 
converges to In other words, for every continuous function (j) with compact support we have 



where dV is the area measure of the horizontal plane Py. In particular, in this example does not converge to 
by even in this weak sense. 

The computation of t„ can be performed as follows. The initial singularity of /"'"(O) reduces to the point 0, 
and hence is a closed subset. It turns out that, for every point x & Py the intensity function by is given by 
by{x) = -{x+ \\x\\v,v) = \\x\\. 

Consider now the domain Dk- The initial singularity of Dk is the set of spacelike lines joining to the 
vertices of the polygon Ck- So when we project on the plane P„, the initial singularity appears as the union of 
k rays si, . . . starting from 0, so that the angle between Sj and Sj+i is li^jk. Let fk'-Py^R the function 
whose graph if the boundary of D/.. On the region Pj of Py that is bounded by Sj and Sj+i, the function 
fk is difFcrcntiablc and the gradient of fk is a unit vector whose angle with Sj is equal to ir/k. The integral 
lines of the gradient are parallel lines that form an angle equal to n/k with both sj and Sj+i. If we denote by 
r{x) the intersection point of the line through x with the singularity, then the intensity function is given by 
bUx) = fkix) - fkirix)) = \\x - r{x)\\. 

Let now Sj be the unique line of this foliation which starts from 0. Clearly, it is the bisector of Pj. If x is on 
the right of s'j then r{x) G Sj. If x is on the left of s'j then r{x) G Sj+i. We now consider the set 

E^ = {xGP,\b'^{x)<\\x\\/2} . 

For each point x & Pj, consider the triangle with vertices at x, r{x) and 0. Note that by{x) is the length of 
the edge joining x to r{x). The sine formula of Euclidean triangles then shows that 

l-vix) = 



sin(7r/fc) 

where </> is the angle at the vertex in the above triangle. Thus, let (j)k be such that sincpk = |sin(7r/A;). Let 
Sj and Sj' be the rays in Pj forming an angle (t)k respectively with Sj and Sj+i. Then, Ej is the union of two 
sectors bounded respectively by Sj and s" and by Sj+i and s"'. 

By the concavity of the function sin in [0, tt/2] we have (pk > 7r/2fc, so for any radius R the area o{ EjnB{Q, R) 
is bigger than ^Area{Pj fl B{0, Rj). Now let us consider the set 

E'' = {x & Py\b';{x) < \\x\\/2} = [jE^ . 

The area of E'' n B{0, R) is the sum of the areas of the E^ n B{0, R), so that 

Area{E'' n B{0, R))>^Y, ^^eo(P,) = ^ 

and, consequently, 

/ {by - b';)dV > [ {by- b^JdV > I . 

Jb{0,R) JB(0,R)nE'' JB{0,R)nE'' ^ 

As Ek is a cone from the origin and the function x ||a;||/2 depends only on the distance from the origin, it 
follows that 



B{0,R) 2 Jb{0,R) 2 

In particular, this shows that by does not converge weakly to by. 
Proposition 4.1. The sequence t*; weakly converges to by/2 in Llg^{Py). 

Proof. Note that since iy{x) < \\x\\, up to passing to a subsequence, there is a weak limit in Ll^^{Py), say J. 
We will prove that J = by/2. This is sufficient to deduce that the whole sequence converges to by/2. 

For this, we consider the following sequence cjk = *duk of 1-forms on Py, where Uk is the function whose 
graph is dDk and * is the Hodge operator. Note that Wk is a L°° 1-form defined in the complement of the 
singularity. 

As duk du at every differentiable point of u and ||(iufe|| < 1, by the Dominated Convergence Theorem we 
have that duk du strongly in Ljg^{Py) as k ^ +oo. This implies that Wfe — )• a; = *du strongly in Llg^{Py). 
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We claim that for any compact supported smooth function / the foUowing formula holds: 
(25) / dfAJuj = ^[ df A . 



Since df Acu — ^dV, this formula implies 



/ ^iij2-J)dV^0, 
J P. dp 



and by a simple density argument we conclude that J 
To prove the claim, first note that 



df AJlj = hm / df A . 

Pv J Pv 

Now, note that on each region Pj of the complement of the singularity, J^ujk is a smooth 1-form, and its 
differential is equal to 

d{iluJk) = dil A Wfc + i^duik 

= di"; A ujk + iyd{^duk) 
= di^ Aujk + i'^AukdV . 

Since on Pj the function Uk is affine, the last term vanishes. Moreover, we have dt^Aujk — di^{gi:a.duk)dV — dV, 
since, on the integral lines of the gradient of Uk, is an affine function with derivative equal to 1. 
Using the fact that tjj A ojk vanishes on the singularity, we finally obtain 



dfAJcj^- fdV, 
Jp, 

which implies Jp df A L^uJk — ^ Jp fdV, and we conclude that 

(26) J df AJuj = hm J df A tjjwfc = - J fdV . 

On the other hand, in oder to compute Jp df A l^ui, note that, on P^ \ {0}, we have the identity 

d{iv Auj) = diy Alu + A doj . 

Now as before diy A uj = dV, but Lydui = LyAudV — ^dV = dV, where the last equality holds since i„ 
this example. So if is the disk centered at with radius e we have 

/ df A LyUJ = - 2fdV - / fiyUJ . 

JPk\B, JPk\B, JdB, 

As cj is bounded, the last term vanishes as e — > so we deduce 



df ALyUj = - / 2fdV. 
Pk JPk 

Equation ([25)1 then follows by comparing this equation with ([25]). □ 

Remark 4.2. Proposition 14.11 makes it clear that the reason why Ly does not converge to l is the fact that Au 
is not concentrated on the singularity, whereas Au^ vanishes outside the singularity. This remark will lead us 
below to introduce the notion of domain of dependence with flat boundary, where this problem is excluded. We 
will then prove fTheorem l4.13p that, for domains of dependence with flat boundary, the convergence does hold. 
We will then show in Section that, without this hypothesis, although the sequence do not converge to Ly, 
it does have a limit, and this limit differs only by bounded factor (attained in the example presented above) 
from by. 

We conclude this section with a simple remark. In the example above we have seen a sequence of domains 
of dependence Z3„ which converges to a domain D, but for which the corresponding sequence of intensities t" 
does not converge to the intensity ly of D. However, in it is clear that 

iy{x) > limsupt"(x) 

n— >-f CO 

This estimate holds in general and is a consequence of two facts: 

• The intensity functions i" of the surfaces of C Dn converge to the intensity function of the C D. 
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• The intensity function of any domain is the infimuin of the intensities of its surfaces of constant 
cosmological time. 
We include a proof of the first statement for the sake of completeness. 

Proposition 4.3. Let he a sequence of domain of dependence converging to a domain D. Denote by 
the intensity function of D and by i„ the intensity function of D. Then for every x d Py we have Lv{x) > 
linisupfc^+^ L^{x). 

Proof. Let us fix e. Denote by the intensity function of the level surface ~ '''k^i^) of th'^ cosmological 
time of Dk- By [J, we know that the sequence of surfaces converges to the level surface = T^^(e) of 
D as k +00. This means that the function : P^, — R, whose graph is H^, converges as fc — )■ 00 to the 
function Ue : P„ — > M which defines H^. By convexity, gradu^(x) converges to gradwe(x). As is given by 



y'l — llgradit^lp 

it follows that t^(x) — i^{x) as fc — )■ +00. Now note that tj(x) < for every k. So passing to the limsup 

we obtain 

lim sup tj; (a;) < i({x) 

and by taking the infimum over e 

limsup 4^(2:) < Lv{x) . 

□ 

4.2. Domains of dependence with flat boundary. Let P be a lightlike plane in R"'^ and denote by g the 
degenerate metric on P induced by the Minkowski metric. We note that P is foliated by lightlike lines which 
are parallel to the kernel of g and denote by P/L be the space of leaves. 

Lemma 4.4. P/L is equipped with a flat metric g that makes it isometric to R" such that g is the pull-back of 
g by the natural projection P — >■ P/L 

Lemma 14.41 implies that there is a natural n — 1 form ujp on P defined as the pull-back of the area form of 
g. This form has the following characterization: if S is any spacelike compact hypersurface in S oriented by a 
future-oriented transverse direction, its area is equal to the integral oiujp on S. 

Now, given a domain of dependence D we consider the 1-form cj on dD defined in the complement of the 
singularity. If p is not on the singularity and P is the unique support plane at x, then uj^ — LOp. 

Regarding dD as the graph of a function u on some fixed spacelike plane Py , it turns out that the identification 
between Py and dD is differentiable at each point where u is difTerentiable, so in the complement of the 
singularities. We can therefore express the from w as a form on P^. 

Lemma 4.5. uj = *du where * is the Hodge star operator of Py. 

Note in particular that w is a L°° form. Moreover, since u is convex, the differential of uj defined as a 
distribution on R" is in fact a positive locally finite Radon measure. More precisely, we have the following 
result. 

Lemma 4.6. du — Au, where Am is a positive Radon measure. 

Proof. Let (u„)„gN be a sequence of smooth convex functions converging to u. Then d * dun — Ait„ — >■ Am 
as distributions. Since the m„ are convex, the Am„ are positive, so their distribution limit Am is a positive 
distribution. Therefore (essentially by the Riesz representation theorem) it is a Radon measure. □ 

Definition 4.7. The boundary of a domain D is called flat, if duj is a measure concentrated on the singularity. 

Let us recall that a Radon measure /i on Py is concentrated on A if /i(Pt, \ A) = 0. It is not difficult to 
check that if the boundary of D is locally given as the union of a finite number of lightlike planes, then it is 
fiat. In Section H75] we will see that a domain of dependence in R^'^ that is the universal cover of a MGHC fiat 
spacetime of genus g > 2 always has flat boundary. The choice of terminology is due to the following lemma. 

Lemma 4.8. Let D be a domain of dependence with a flat boundary. Suppose that A is an open subset of dD 
which does not meet the singularity. Then A is contained in a lightlike hyperplane. 

Proof. As Iless(M) is a measure with values in the positive definite quadratic forms, one finds that if Am — 0, 
then Hess(u) is zero on A. This implies that —u is a convex function too, so u is an affine function, i. e. the 
graph of u\a is a plane. Since it is contained in the boundary of D, this plane must be lightlike. □ 
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Remark 4.9. If is a domain of dependence obtained as the intersection of the future of a locally finite family 
of lightlike planes in M"'^, then its boundary is clearly flat according to the previous definition. Moreover, 
Lemma 14.81 shows that if D is a domain with fiat boundary and the initial singularity is closed in Z?, then D is 
the intersection of the future of a locally finite family of lightlike planes. 

However, in the next subsection we will show some interesting examples of domains of dependence with 
flat boundary which are not of polyhedral type. This depends on the fact that in those examples the initial 
singularity is not a closed subset. In fact, in many cases the singularity is dense. 

If is a domain of dependence with flat boundary, the following proposition holds. This proposition will be 
the key ingredient in the proof of the stability of the intensity functions. 

Proposition 4.10. Let v be a timelike unit vector. If D is a domain of dependence with fiat boundary then for 
every compactly supported smooth function f on Py the following identity holds: 

df A M ^- fdV . 

The proof of the proposition is based on the following lemma, which is also valid if D does not have flat 
boundary. 

Lemma 4.11. Let : P„ — >■ M 6e the function whose graph is the level surface and let be the intensity 
function of . Then at every point where gradwe is differentiable we have 

(27) < (grad/,e,gradu,) < 1 , 

LfA{ue) < n — llgradwelp . 

Proof. Since the vector gradu^ + w is orthogonal to the surface H^, the unit normal future-oriented vector is 
obtained by normalizing it as 

:(graduc+w). 



y/1 - llgradu^ll^ 



which implies 



\/l - ||gradwe||2 
Thus at every point where ||gradMe|| is differentiable we have 

(gradte,gradu,) = _ [[gradu p)3/2 (Hess('Ue)grad'Ue, grad^e) , 

which is positive by the convexity of u^. 

To prove the estimate from above we use a comparison argument. Consider the retraction of the point 
q = X + Ue{x)v on the singularity, say r = r{q). Note that on /"""(r) n D the function f{p) = r,p — r) 

is less than than the cosmological time. It follows that the level surface H' = /^^(e) is contained in the closure 
of the future of H^. Moreover, since T{q) — f{q) = e those surfaces are tangent at the point q. 

Note that H' is the graph of the function 



h: Py^ M, h{x) =c + v/||a;-r||2 -f 

where c G M and f d Py are determined by the orthogonal decomposition r = f + cv . 

From the fact that H' is tangent to iZ at g and that it is contained in its epigraph, one deduces 

• u^{q) = h{q), 

• gradu<;(g) = grad/i(q), 

• HessMe(g) < IIess/i(g), 

where the last inequality implies that the difference is a positive definite matrix. In particular, we find that at 
the point q 

(gradte,grad'Ue) < 7- ^ ^ — (Hess(ft)grad/i, grad/i) . 

(1 — llgraduclp)'^/'^ 

Now an explicit computation shows that 

grad/i = {x — r) , Hess/i — {Id — grad/i ® grad/i), 

h — c h — c 
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which inipUes that, stih at the point q, the fohowing inequahties hold 

(grad(t,),grad(u,)) < \u2\3/2tu T(llgi'ad(/i)|p - ||grad(/i)f ) 

(1 — ||grad('Ue)|pj''/"'(/i — c) 

< 7; n '^Mi2A3/2^ r(llgradK)f - ||grad(we)||'' 

(1 — ||grad('Ue)|p)-^/^('u — c) 

Using ||gradue|| < 1 and the identities 

Ue(g) - c = -(q - r{q),v) = L^{q) = 



^1 - ||grad(ue)|l 
we obtain 

(gradte, gradUf) < |jgraduc||^ < 1. 
To prove ([28|. it is then sufheient to note that 

Aue{q) <Ah = -^{n - llgraduein • 



□ 



Proof of Proposition [7^T0[ Let be the intensity of the surface H^, let : M be the C^^^-function whose 

graph is the surface H^, and set = *dUf. Then — > in Lj^^. Moreover ui^ is a Lipschitz form. Analogously, 
we find that is a Lipschitz function since 

VI - llgradweP 

It follows that for every smooth function with compact support / we have 

(28) / df A {l,uj,) = - / fdi, Auj,- [ fi,Au,dV . 

As te \ /, pointwise and / has compact support, there exists a constant C such that \fie\ < C for e < 1. It 
then follows by the Dominated Convergence Theorem that 

fi^Au j fiAu = , 

where the last equality holds because t is zero on the singularity and Au is concentrated there. 
On the other hand, we have Au^dV AudV as measures, which implies 



j fi^u 


< C 


[ Au,dV~ [ Au 









Jk Jk 





As t vanishes on the singularity, it follows that the last term on the right hand side of ([28]) converges to 0. To 
conclude, it is sufficient to show that 



fdbf Aujc ^ fdV . 
'P„ JPv 

Now the 2-form dtf A uj^ is equal to dL^{gTSLdug)dV . So it is sufficient to prove that — (gradte, gradug) weakly 
converges to 1 in Lj^^. 

First we consider the case where D is the intersection of the futures of a finite number of lightlike hyperplanes - 
we will describe this situation by saying that D is "finite" . Then the singularity is a finite tree and, in particular, 
it is closed. In this case, given a point x G we denote by ro(x) the starting point of the lightlike ray through 
q{x). Then the restriction of the map Tq on each region E of dD \ T is a smooth projection and satisfies 

Lv{x) = -{q{x) - ro{x),v) . 

A simple computation shows that gradtu = gradit on E, which implies 

(29) (gradt„ , gradu) = 1 . 

On the other hand, if ro{x) lies on the interior of a segment e of T, it is not difficult to check that gradi^ (x) -> 
gradt(a;). Indeed consider the domain D defined as the future of the spacelike line which contains the segment 
e. Note that D D D, and e is contained in the singularity of D. Thus if f : _D — > e denotes the retraction on 
the singularity, then that f~'^(e) = r~^(e) = U, and the cosmological time of D coincides with the cosmological 
time of D on U. Then the intensity iy and around q can be computed by considering the domain D instead 
of the domain D. In that case an explicit computation shows that the convergence to (|29)) holds. 
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In particular, the function (je converges to 1 almost everywhere. Since Lemma 14.111 shows that is bounded 
by 1, the Dominated Convergence Theorem implies that converges strongly to 1 in Lj^^. So the proposition 
is valid, whenever D is finite. 

Consider now the general case. By Lemma 14.111 we have < 1 at every point. So for every sequence — >■ 0, 
up to passing to a subsequence, we can take the weak limit in Lj^^. That is, there exists a function g such that 

fgdV= liin / fg.^dV , 

n— >-+oo J 

and ||.g||L°° ^ 1- Note in particular that for every smooth function with compact support we obtain 

dfALLO^-f fgdV 



Clearly, the same formula holds for any function / which is the limit in the Sobolev space W^'^{K) of smooth 
functions with compact support. 

To prove that g = 1 we use an approximation argument. Let Dk be a sequence of finite domains of dependence 
converging to D, and let : — R be the functions whose graph is dDk- Consider a smooth function 
(j) : [0, +00) — )• [0, +00) which is decreasing and such that its support is contained in [0, M], and define 

fk{x) = (l){uk{x)) , f{x) = (l){u{x)) . 

Then the functions fk are with compact support and are limits of smooth functions with compact support. 
They satisfy 

dfk = (l>'iuk)duk , df = (j)'{u)du . 
As holds for finite domains we obtain for every k 

-(j>'{uk)ikdV = [ cl>{uk)dV , 



and with the inequality (j)' < we deduce 

\^'{uk)\ikdV ^ I q>{uk)dV 



As Ufc — ^ u uniformly on compact subsets and 4>{uk) is zero outside a compact region, which is independent of 
fc, we obtain 



(30) j \c^' {uk)\ikdV I ^iu)dV 
On the other hand, Fatou's Lemma implies 

(31) lim / \(j}'{uk)\LkdV < I limsup(|0'(Mfc)|ife)dy 



Note that by Proposition 14.31 the right hand side is less than J\(f)'(u)\LdV which is equal to / 4>{u)gdV. So 
comparing ([30]) and (|3T|) yields 

' \cl)'{u)\dV < j W{u)\gdV, 

and since g{x) < 1 almost everywhere, it follows that g{x) = 1 almost everywhere. □ 

Remark 4.12. In other words, Proposition 14.101 states that LyUj is a primitive of the volume form of Py in a 
distributional sense. 

4.3. Stability of the intensity for domains with flat boundary. We can state now the main stability 
theorem of this section. 

Theorem 4.13. Let Dk be a sequence of regular domains converging to D. If the boundaries of Dk and D are 
all flat then — >■ strongly in Lj^^{Py). 

Proof. Consider the convex functions Ufe : — )■ R whose graph is identified with dDk ■ Note that the sequence 
Uk converges to u uniformly on compact subsets. Moreover, since duk converges to du almost everywhere and 
||c?u/c|| < 1, it follows from the Dominated Convergence Theorem that duk — t- du in L^{Py), which implies 
Wfe ^ w in L^{Pv). 
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We may assume that Uk > for every k. By this assumption, < Uk{x) for every x. So the functions 

are uniformly bounded on compact subsets. Consequently, there exists a weak limit J of the sequence ijj, that 
is a function J with 

(32) / fitdV ^ / fJdV 



for every compactly supported continuous function /. 

Now we claim that for every compactly supported smooth function / we have 

(33) J df/\ ( Jw) ^- J fdV = J df A M . 

From this claim it then follows that Jlj — l^lu almost everywhere and hence J = by. 
To prove the claim, first note that the left-hand side of (1551) can be rewritten as 

(34) f dfA{jLu)^f (J-4-)d/A^+ / L^df A[uj-LOk)+ I dfAiiluk). 

J J J Pjj J Pjj 

By Proposition 14. 10[ the last term is equal to the right-hand side of ([55)1 and hence independent of k. 

As uJk converges to uj in L]^^{Py), the second term in ([M)) vanishes as — >■ -l-oo. (Note that the functions 
L^df are uniformly bounded in L'^^). 

Finally by the density of continuous functions in Lf^^, J is the weak limit of in Lf^^, that is, <\'6'2\i holds 
for every L^-function defined on some open subset with compact closure. In particular, the first term on the 
right in p4p also vanishes as fc — > +oo (indeed df Alu — gdV for some compactly supported L°°-function g). So 
letting k go to -t-oo in we deduce that (I55| holds, and the claim is proved. 

To show that is a strong limit of tjj we use the following remark: 

Vx e dD : Ly{x) > limsup tj;(a;) . 

This holds due to the following: if is the maximal lightlike ray in dDk that contains the point x + Uk{x)v, 
then the sequence converges to a lightlike ray r contained in dD^ which contains the point x + u{x)v. So the 
same argument as in the proof of Proposition 14.31 may be used. 

Let us now fix an open subset A C Pi, with compact closure. Then there is a constant M > such that 
and tjj are bounded by M on A. This implies that the functions — (tjj)^ are uniformly bounded and allows one 
to apply Fatou's Lemma. It follows that 



/ -{ivf < I -limsup(tS)2 < hminf f- / {ilf 

J A J A \ J A 



which implies 

\Vv\\l^{A) > IhTlSUp ||t^||L2(^) . 

as the sequence converges weakly to t„ in L^(A), this estimate implies the strong convergence in L'^{A). The 
L^-convergence on a compact set then implies the L^-convergence. □ 



4.4. Uniform bounds for domains with non-flat boundaries. We have seen in Section H73l that if D is 
a domain of dependence with flat boundary, and {Dk)k&i is a sequence of domains of dependence with flat 
boundaries converging to D, then the intensity of the Dk converges to the intensity of D. Suppose now that D 
is any domain of dependence (not necessarily with flat boundary) and that {Dk)k&i is a sequence of domains 
with flat boundaries converging to D. We know (Proposition 231) that the intensity of D is at least the limsup 
of the intensities of the domains but the example in Section [4. II shows that equality does not always hold. 
We will now see that the opposite inequality does holds, albeit with a multiplicative constant. 

This result is important in view of the numerical computations in the following sections, which allow one 
to visualize the intensity seen by an observer in an MGHFC manifold of dimension 3-1-1. These computations 
are done by approximating the corresponding domain of dependence by a sequence of finite domains with flat 
boundaries. However it might happen that the limit domain does not have flat boundary. The computed 
intensity function then coincides with the limit of the intensities of the finite domains (with flat boundary) , and 
can differ from the actual intensity of the limit domain. However, Theorem 14.141 then ensures that the actual 
intensity is at least equal to the computed limit intensity, and that it is at most three times this computed 
intensity. 

Theorem 4.14. Let D C R"'^ he a domain of dependence, and let {Dk)k&i ^2 a sequence of domains of 
dependence with fiat boundaries converging to D. Let t„ he the intensity of the boundary of D with respect to a 
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unit timelike direction v, considered as a function on Py, and let iJ^ he the intensity of Dk- Then the sequence 
converges in Lj^^ to a limit bum, and 

Him — by ^ niii^i . 

In Section 231 we will show that all domains D that arise as universal covers of 2+1-dimensional a MGHFC 
spacetimes have flat boundaries. Consequently, in that case the inequalities in statement of the theorem can be 
improved to an equality. In higher dimensions, it appears unlikely that the boundary of the universal cover of 
a GHMC fiat manifold always has a flat boundary. It is quite conceivable that for such domains, the inequality 
can be improved to 

blim — by (71 ^)blim : 

but we do not pursue this question further here. 

Proof. We proceed as in the proof of Theorem 14.131 and only indicate the steps that differ from that proof. As 
already noted. Lemma [4.111 still holds, but differences occur in the proofs of Proposition 14. 101 and of Theorem 
14.131 As we want to obtain inequalities on by we consider a positive test function /. Equation (pS)) still holds. 
However, the inequality is weakened to 

0< J fbyAu<{n-l) J fdV . 

The last inequality descends by the estimate when one takes the limit e — > and uses that Auc Au as 
measure and that ie \, i^. Following the proof of Proposition 14. lOl we then obtain the inequality 

[ fdV < [ dfA M < - [ fdV . 
In the proof of Theorem 14.131 Equation ([55)1 is therefore replaced by 

[dfA {Ju) - - / fdV , -n [ fdV< [ dfA (byuj) <~ [ fdV . 
J p^ J p^ J p^ J p^ J p^ 

The rest of the proof of Theorem 14. 131 goes through and leads to the statement. □ 

4.5. Universal coverings of MGHFC spacetimes in dimension 2+1. Let M be a 2+1-dimensional MGHC 
flat spacetime of genus g > 2. We know that the universal covering of M is a domain of dependence D C M'^'^. 
We will assume in the following that M is not a Fuchsian spacetime. This means that D is not the future of 
a point, or, equivalently, that the holonomy representation of M is not conjugate to a linear representation in 
SO {2, 1). In this subsection we will prove that the boundary of D is flat. This is a consequence of the following 
geometric property of the boundary of D. 

Proposition 4.15. Identify the boundary dlP with the set of lightlike directions in M.^'^ and let D* be the 

subset of dip consisting of lightlike directions parallel to lightlike rays contained in dD. Then the set D* has 
Lebesgue measure zero in dM? . 

We know by the work of Mess [TO] that the linear part of the holonomy representation of M defines a Fuchsian 
group r, which determines a hyperbolic surface S — H^/F. Moreover, there is a measured geodesic lamination 
A on such that M is obtained by a Lorentzian grafting on the Minkowski cone of S. Denote by A the lifting 
of A to the universal covering H^. We say that a point ^ e ffM? is nested for the lamination A if for some point 
V £ H^, the intersection of the ray joining w to ^ with A is +00. 

Lemma 4.16. //^ is a nested point for A then 

• the point ^ is not the end-point of any leaf of X, 

• the intersection of any ray ending at with A is +00. 

Proof. Let us consider the upper half-plane model of H^. Without loss of generality we may assume that ~ 00. 
Suppose there is a leaf I of A ending at ^, and take any compact ray tq joining a point v e H'^ to I. Now any 
sub-arc of the ray [f,^) can be deformed through a family of transverse arcs to a subarc of tq. This implies 
that the intersection of any subarc of [v,^) with A is uniformly bounded by the intersection of rg with A. This 
proves that <^ is not nested. 

For the second statement, consider a point vq e H-^ such that the intersection of [wqjC) with A is +cx). Take 
a family of leaves In meeting [vq.^) at a point Vn — >■ With the first statement, it is easy to check that In 
bounds a neighborhood [/„ of ^, and that {?/„} is a fundamental family of neighborhoods of ^. 
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In particular there is a leaf, say li, cutting both [wqiC) a point vi and [w,^) at a point wi. Then every 
leaf of A cutting [vi,£_) must cut also [wi,^) This implies that the intersection of [w,^) with A is bigger than 
the intersection of [wi,^) with A, which is clearly infinite. □ 

We will see that the set D* does not contain any nested points, so the proof of the proposition then follows 
from the following fact of hyperbolic geometry 

Lemma 4.17. Almost all points in diP are nested for A. 

Proof. The proof is based on the ergodicity property of the geodesic flow on S. For {x, v) in the unit tangent 
bundle of S let r(x, v) the geodesic ray {exp^{tv)\T > 0}. Consider now the following subset of T^{S): 

Bn = {{x, v) € T\S)\i{r{x, v), \) < n} . 

We claim that Bn is a set of measure zero for the Liouville measure. 

Before proving the claim, let us show how the claim proves the statement. Indeed we get that the measure 
of the set B = [_\ Bn is zero. Let B C T^H^ be the lifting of B on the universal covering. By definition we have 
that B is made of pairs {x, v) such that the endpoint of the ray exTp^{tv) is not nested. In particular, if E is the 
complement in SH^ of nested points, the Liouville measure of B can be computed as as 

JK 

where K is a fundamental region and fix is the visual measure from x. As the measure of B is zero, it immediately 
follows that E is a, set of measure zero. 

It remains to prove the claim. Let <j)t denote the geodesic flow on T^S. Clearly we have 

MBn) C Bn. 

More precisely, t < s implies 0f(i?„) C (f)s{Bn). It follows that i?„ — Ut>o'^t(^) ~ UfceN'^*:(^n) ^ subset 
invariant by the geodesic flow. Moreover its Liouville measure is equal to 

lJL{Bn) = inf ^(0fe(B„)) = fl{Bn) 
k 

where the last equality holds because fi is invariant by the geodesic flow. 

By the ergodicity of the flow, we have either fi{Bn) = or fj,{T^S \ Bn) = 0. In order to prove that the latter 
is not true, it is sufficient to prove that the complement of i?„ contains a non-empty open subset. 

First note that if (x, v) corresponds to a closed geodesic which intersects A, then the intersection of A with 
the ray exp^(tu) is +oo. In particular {x,v) ^ B„. 

Moving X on the ray, we may assume that it is not on A. Now take Af > so that the intersection of A with 
the segment r = {exp{tv)\t e [0, M]} is bigger than 2n and exp^(Mw) is not on the lamination. 

We want to show that a neighborhood of {x,v) is contained in the complement of Bn. Indeed if {xk,Vk) 
converges to {x,v), then the intersection of A with the segment rk — {exp^^{tvk)\t S [0, M]} converges to the 
intersection of A with r. So for k sufficiently large, {xk,Vk) does not lie on Bn. D 

In order to relate nested points with points in D* we need this technical lemma of Lorentzian geometry. 

Lemma 4.18. If R is a lightlike ray contained in dD which is maximal with respect to the inclusion, then there 
is a sequence of points r„ on the singularity T which converges to a point on R, and a sequence of .spacelike 
support planes P„ at r„ which converges to the lightlike plane containing R. 

Proof. Let v be any future oriented timelike vector. Let g be a point of r and consider the segment of points 
Qe — q + ev for e G [0, 1] and the path on E given by r^ — r{q^). 

Note that this path is contained in the closure of D n I^{qi), which is a compact region of R^'^. Thus there 
exists a sequence e„ — >■ such that r^^ converges to some point f. 

We claim that f is contained in R. In order to prove the claim, note that sequence of vectors q^^ — r^^ 
converges to g — f . Since they are timelike, their limit cannot be spacelike. But dD being achronal forces q — f 
to be lightlike and the segment [f, q] to be contained in dD. As the lightlike plane P containing i? is a support 
plane for _D, it follows that [f, g] is contained in this plane, so in particular is on R. 

To construct the sequence of lightlike support planes, it is sufficient to set P„ to be the plane orthogonal to 
q^^ — r^^^ passing through r^^ . □ 



We are now ready to prove Proposition 14. 15) 
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Proof of Proposition \4-15\ We will prove that if ^ G D* then ^ is not nested. Assume by contradiction that ^ is 
nested, and let i? be a ray parallel to ^. By Lemma [4. 181 there exists a sequence of points r„ on the singularity, 
converging to a point on R and a sequence of spacelike support planes P„ converging to the lightlike support 
plane containing R. Let Un be the unit timelike vector orthogonal to P„. Clearly we have that m„ — ^ ^ in H-^. 
By Mess' construction [TH], we have 

rn-ro^ Wn{x)d^jLx, 

where c„ is the segment joining uq to u„ and Wn{x) £ M^'^ is defined to be if x is not in the support of A and 
is the unit tangent vector at x orthogonal to the leaf through x and pointing towards u„ otherwise. 

As by the hypothesis ^ is nested, the ray r joining uq to ^ transversely meets the lamination. In particular, 
by changing uq to a point on r H A we may assume that uo is on the lamination. Let e be the unit vector at mq 
orthogonal to the leaf through uq and pointing towards ^. We claim that if ^ is nested then 

(r„ - ro,e) +co, 

which contradicts the assumption that the sequence r„ converges in Minkowski space. 

First note that since Un — >■ ^, we may assume that m„ is on the half-plane bounded by Iq and containing ^. 

Thus if Hs a leaf that intersects c(Mo,Mn), I disconnects Iq from ^, and the scalar product of vectors e and 
w{x) is positive. Since the corresponding geodesies are disjoint, the reverse of Schwarz inequality holds, that is, 
(w(x),e) > 1. This implies 

(r„ - ro,e) > i{\,Cn) ■ 

Let us prove that t(A,c„) — > +oo. The reason is that for every leaf I of A cutting the ray [uq,^), Un is 
definitively contained in the region bounded by I containing ^. So for every point x on the segment [uo,0; foi' 
n sufficiently large, we have i{i>,Cn) > [uq,x]). 

Since we are assuming that J^([wo:C)) ^ we can choose x so that [uo,x]) is arbitrarily big, so the 

conclusion follows. □ 

Let us fix a unit timelike vector v, and let u : P„ — >■ _D be the convex function whose graph is the boundary 
of D. Note that if e is a unit vector in P^,, then e + w is a lightlike vector. In this way, the unit circle in Py 
is identified to dH^ by the map e t-^ [e + v]. Under this identification, the subset D* corresponds to the image 
of the map 

6 : Py\T ^ S\ 6{x) = grad(u)(a;) . 
Fix a unit vector e in Py and take linear orthogonal coordinates {x,y) on Py such that dx = e, and consider 
the restriction of the function u on each line parallel to e. That is, for e M, let Uj^ : R — > M be defined as 
Uy{x) = u{x,y). 

Note that whenever [x, y) does not correspond to a point on T, then Uy is differentiable at x and 

{uyYix) = (gradu(a;,y),e) 

By Proposition 14. 151 at those points, the derivative takes value in a subset of zero measure of M. 

Now to prove that the boundary of D is flat we will proceed in three steps. 
Step 1. We will prove that for a generic choice of the vector e, if Uy is differentiable at x, then {x,y) does not 

correspond to a point on T. In particular, the derivative {uy)' takes value in a subset of zero measure 

of R. 

Step 2. We will use this fact to show that for every y, the measure (uy)" is atomic with support on Tn Ry. 
Step 3. Using a disintegration formula for dxxU in terms of the family of measures {uy)" we conclude that this 
measure dxxU is concentrated on T. 

Lemma 4.19. There is a subset A of such that: 

• If e £ A then, for every y g R, the points x where Uy is differentiable are exactly the points such that u 
is differentiable at (x,y). Moreover, at those points, 

(uy)' = {gradu{x,y),e) . 

• The measure of \ A is zero. 

Proof. Let A be the set formed by vectors e such that the geodesic in starting from v with direction e does 
not meet any leaf of A orthogonally. We will prove that A fulfills the requirements of the statement. 

First let us prove that the only differentiable points of Uy correspond to differentiable points of u. By 
contradiction suppose that Uy is differentiable at a point x so that {x, y) corresponds to a point on the singularity. 
Up to translation we may suppose that x ~ y ~ u{x, y) ~ 0. 
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As the point is on the singularity, there are two hghthke planes Pi, P2 through 0, which are support planes 
for D. Let V be the vertical plane containing e and v. Note that Pi n F and P2 n 1^ are support lines for dD n V 
at the point 0. As dDr\V corresponds to the graph of Uy, and we are assuming that Uy is differentiable at a; = 0, 
those support lines must coincide, Pi n = P2 n This implies that V must contain the line r = Pi n P2. 
Note, however, that this line is spacelike, and its dual geodesic in is a leaf of / G A. On the other hand 
V n is the geodesic g starting from v with direction e, so the condition implies that g meets orthogonally I, 
contradicting the choice of e. 

It remains to show that the complement of A in is a set of measure zero. Note that if e G A then — e is 
also in A, so we may regard ^ as a subset of the projective line P{Pv). 

We will argue as follows. For any geodesic I of H^, let e{l) be the unit tangent vector at v such that the 
geodesic exp„(te(^)) hits orthogonally I. Note that e{l) is defined up to the sign, so it should be considered more 
properly as an element of P(P„). The complement of A is the set of unit vectors e{l) where I is a leaf of the 
lamination A. 

Let us now fix any ray r starting from v and define 

E,. = {e{l) \ I is a leaf of A hitting r } . 
Note that if Vn is a dense subset of the ray from v we clearly have 

\jEr„= P{P,)\A . 
n 

So in order to argue that the measure of the complement of A is zero, it is sufficient to show that Er has measure 
zero. 

Now on r n A we may define a vector field w by taking for 'w{x) to be the unique vector orthogonal to the 
leaf I through x. Note that e(Z) coincides with the orthogonal projection of w{x) on Py up to renormalization. 
By a classical result [TT], the field w can then be extended to a Lipschitz vector field, still denoted w, on the 
whole line. In particular, we obtain a map 

e : r ^ P(P,) 

by defining e{x) to be the projective class of w{x). It is not difficult to show that this map is locally Lipschitz 
and, by definition, = e(r n A). As the measure of r n A is zero, this concludes the proof. □ 

Lemma [4.191 concludes the proof of step 1. In particular, note that if e is on the set A, then for every y eM., 
the derivative function Uy takes value on the set 

{(gradu(a::,y),e)} 

which by Proposition 14. 15] has measure zero. The proof of step 2 is then based on the following simple lemma 
on convex functions. 

Lemma 4.20. Let u : K — >■ R &e a convex function. Suppose that the measure of the image u' : R — > M is zero. 
Then u" is an atomic measure, and its support coincides with the set of discontinuity of u' . 

Proof. By the standard theory of convex functions the measure u" can be split as the sum of a measure fx 
without atoms and an atomic part, say v, with v = cikSq^, where Sq^ is the Dirac measure concentrated on 
qk and ^ Ok is an absolutely convergent series. 

Now we claim that the measure of the image of u'{[a,b]) is equal to /i([a, 6]). Indeed, note that u is not 
differentiable exactly on the points {qn\n G N}. Moreover at every point there exists the left derivative and the 
right derivative that can be expressed as follows. Assume that is a differentiable point of u, then 

u[{x) ^ u' (0) + (j){x) + On , u'^{x) = u'{0) + (t){x) + ^ a„, 

where we put (f>{x) — /i([0, x]). 

If u is differentiable at some point x, then the values that u' takes on the interval [0, x] can be described as 

u\[0,x]) = [u\0),u'{x)]\ y /„ 

where /„ — [u'i(qn), u'r{qn)] is the interval between the left and right derivative at q„, which are pairwise disjoint. 
So the measure of this set it'([0,x]) is given by 

(j){x) ^ a„ - ^ a„ = (j){x) . 

g,ie[0,a:] g„e[0,a:] 

By the assumption, (j){x) = for any x, so /i = 0. □ 
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Finally, in order to prove step 3 we need the following disintegration result of the measure d'^^u in terms of 
the measure [uy)" . 

Lemma 4.21. Let x, y be coordinates on and consider a convex function u. For every y g denote by 
Uy the convex function x i— > u{x,y). If d^xU denotes the second derivative of u along the x axis (which is a 
Radon measure on M^j and (uy)" denotes the second derivative of Uy (which is a Radon measure on then 
the following formula holds 

/ f{x,y)dxxu= dy f{x,y){uy)" 

for every bounded Borelian function f with compact support. 

The proof of this analytical Lemma can be found in 21 (Theorem 1.3 formula (1.31)) for the wider class of 
bounded Hessian functions. We are ready now to prove that the boundary of D is flat. 

Proposition 4.22. // D is the universal covering of a MGHC flat spacetime of dimension 2 + 1, then its 
boundary is flat. 

Proof. We will prove that Hess(u), considered as a matrix- valued measure on P^, is supported on T. 

Indeed it is sufficient to prove that there are three independent directions 61,62,63 such that D^, g. (u) is a 
measure supported on T. As the subset A from Lemma [4.191 is dense, it is sufficient to prove that ^(u) is 
supported on T for 6 € A. 

If X, y are the Cartesian coordinates on Py such that e — dx, then g(u) coincides with 9^^(u). So we have 
to prove that if / is a measurable function which is zero on T, then 

J fdlM = . 

We may compute the integral above using Lemma 14.211 which yields 

/ /C- / dyj{fy){uy)" . 

Now, by Lemma [4.201 and Lemma [4.191 [uy)" is supported on the discontinuity of u'y which corresponds to 
points X such that {x,y) in on the projection of the singularity. It follows that fy is zero on the support of 
{uy)" , and hence the integral is zero. □ 

5. The intensity function for SPACETIMES constructed from measured geodesic LAMINATIONS 

The aim of this section is to understand to what extent and under what conditions an observer in a 2 + 1- 
dimensional domain of dependence can reconstruct the geometry and topology of the ambient space from his 
observation — either at one time or over a fixed time interval — of the intensity of the signal emitted by the 
initial singularity. In particular, we investigate this question for a domain of dependence, which is the universal 
cover of a (non-Fuchsian) MGHCF spacetime M containing a closed Cauchy surface of genus g. 

A basic remark, which somewhat complicates the statements and the analysis below, is that the observer can 
only "see" the universal cover of Af , so he can in no way distinguish M from any of its finite covers. In other 
words, the observer can only determine the largest discrete subgroup of Isom(R^'^) compatible with the signal 
she observes. Moreover, he can only be certain to have determined correctly the fundamental group of M if he 
knows the genus of S, since otherwise it remains possible that his spacetime is topologically a finite cover of M, 
with a flat metric which is "almost" lifted from a flat metric on M , with only a small change in a region not 
visible by her. 

In Section 15.11 we study the relationship between the intensity function measured by an observer and her 
cosmological time (see Proposition 15. 6p . Then in Section 15.21 we show (see Proposition I5.13P that an observer 
in the universal cover of a non-Fuchsian MGHCF Minkowski spacetime can reconstruct in flnite eigentime the 
geometry and topology of the space, if the genus of the Cauchy surface is known to him. In Section 15.31 we 
briefly explain how those arguments can be adapted to higher dimensions. 

5.1. Estimating the cosmological time from the intensity function. We consider a domain D and an 
observer in D given by a point p E D and a future directed timelike unit vector v G H^. We consider the 
associated rescaled intensity function 

P^ u : -S*^ Mq 

defined as in the previous section. We will be mostly interested in the case where D is the universal cover of 
a spacetime constructed by grafting along a measured geodesic lamination. In this case, the observer will see 
a division of the circle into intervals, on which the intensity function behaves like that of a spacelike line, and 
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intervals in which the intensity function behaves hke the one of a hght cone. The former correspond to the edges 
of the singular tree of D, the latter to its vertices. It follows from the results in the previous sections that the 
intensity function is analytic on the segments of the circle that correspond to the edges, while it is generally not 
analytic on the segments that correspond to the vertices. In general, the segments of the circle that correspond 
to vertices of the singular tree form a Cantor set. 
We define the maximum intensity as 

p^f-" = sup p^Ju) 
To understand its properties, we consider again our two main examples. 

Example 5.1. Consider a cone D — I^{q) and an observer characterized hy p E (q) , v e H^. Then 

D,rnax 

Pp,v 2 ' 

where T{p) ~ \p — q\ is the cosmological time of the observer and S = d^2 {v, grad^ T) is the hyperbolic distance 
between v and grad^T. This follows from Equation (|4]). 

Example 5.2. Consider the future D = of a spacelike line I and an observer with p G v G H^. It 

follows from Equation (j20l) that 

where T{p) is the cosmological time of p and S, ^ are defined as follows. Denote by I the geodesic in M that is 
stabilized by the PSL(2,M.) element that fixes the direction of I. Then gradp T defined a point on I, 5 is the 
hyperbolic distance of v and I and ^ the distance of the projection of w on ^ from the point in defined by 
gradp T. 

We will now determine an estimate on the cosmological time. The central ingredient is the following propo- 
sition. 

Proposition 5.3. The maximum intensity is given by 

^D.jnax _ g^pj^^^ y — x) : X , y E (p) H D , y — X future directed and lightlike} . 
Proof. By definition, we have 

pD-jnax ^ sup{^^^ y — x) : X , y € (p) (1 D , y — X future directed and lightlikc} . 
To show the opposite inequality, we choose sequences (a;„)„gN, (yn)neN with x„,y„ e I {p) ^ D that satisfy 

(35) lim {v, yn — Xn) = supKw, y — x) : x,y E {p) C\ D,y — x future directed and lightlike} . 

As the intersection I^{p) H Z) is compact, there exist convergent subsequences (a;„j.)fegN, (l/nfc)feeN with limits 
e / (p) n D, yn^^ -?■ y & I {p)r\ D. 
If the segment [x, y] is extensible, i. e. if there exist x', y' G I {p) H D with [x, y\ C [x\ y'], then we obtain a 
contradiction to ((55)) . Therefore [i, y] is inextensible, x lies on the tree and y G dD n I^{p)- This implies that 
there exists a 9 E such that the past directed lightlike ray starting at p that is defined by intersects dD in 
y and {v,y-x) ^ p^,,{9). □ 

An immediate consequence is that if an observer moves along a timelike geodesic then the maximum intensity 
is increasing with time. More generally. Proposition 15.31 allows us to give estimates for the maximal intensity 
of domains that are contained in each other. 

Corollary 5.4. Let D, D' be domains with p e D C D' . Then for all w G H^.- 

D,max ^ D\max 
rp^v — rp.v 

In particular, we can estimate the maximum intensity function for any domain. 

Corollary 5.5. For a domain D and any observer characterized by p E D and v G H^, the following inequalities 
hold: 

(36) p^r"" > sup , p^r"" < , inf . 

Q^J- (r))riD ^ spacetike tine 

DCI+(l) 

By applying this corollary to a general domain and using the results of Examples 15.11 and 15.21 we obtain the 
following statement. 
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Proposition 5.6. Let D be a domain with an observer characterized by p E D and v G H'^. Then the following 
inequalities relate the maximum intensity and the cosmological time: 

(37) ^ < p^;"''" < T{p) coshdH(^^, grad^ T) . 

Proof. From the first inequahty in CoroUarv 15.51 and Example 15.11 we have 

P?:r^> sup p^;;'''"'^^-^ sup Ip-gle'^-C".^) > i sup \p-q\. 
q£i-(p)nD ^ qei-{p)nD ^ qei-{p)nD 

By definition of the cosmological time, T{p) = sup^gj-^p-j^^ |p — gj, which proves the first inequality. 

To prove the second inequality, we use the fact that for any lightlike vector ^ and any two future-directed 
timelike unit vectors m, u, we have 

(38) \{v,0\<2\{u,v)\\{u^0\- 

This can be seen as follows: after applying suitable elements of 50(2,1)+, we can suppose that the vectors 
£,,u,v are given by 

(1 \ / cosh a \ / " \ 

v=i sinha ^= 6 with = b"^ + . 

/ V y \c J 

This yields 

|(u,^)| = |a|, Kw,^)! = jacosha — bsinhaj < 2|a| cosha , |(u,ti)| = cosha , 
and proves By combining with Proposition 15.31 we obtain for all p € D and u,v eM'^ 

For u = gradp T this yields 

p^;— < 2coshdH(i.,grad,(r))p^g7,7^^ . 
For the future of a spacelike line, we have from Example 15.21 

/+((), max _ T{p) 
^p.gradp T — 2 ' 

because the parameters 5,^ in Example 15.21 vanish . This implies together with Corollary 15.51 that 

p^r"^ < 2 coshdH(«, grad,(r)) inf p^^'lrr ■ 

I spacchkc line ^ ^ p 

By definition of the domain D there exists a point q in the tree with T{p) = \p — q\ and two lightlike support 
planes that contain q. Let I be the line obtained by intersecting these support planes. Then the cosmological 
time T{p) of p with respect to I'^il) and its gradient grad^ T at p coincides with the cosmological time T{p) 
with respect to D and its gradient gradp T. This implies 

_T{p) 

, 1- ^P,gi-ad„T ~ Pp,sra.d (T) ~ o 

i spacolikc Imc '^"^ i' f:B^'^^p\J- / / 

DCl+il) 

and proves the claim. □ 

5.2. Reconstructing the holonomy from the intensity function. From the intensity function on the 
circular segments, the observer can reconstruct the relevant data (position of the edges and vertices, his geodesic 
distance from the edge segments and vertices), but only for the pieces of the singular tree he sees. If the observer 
is very close to the singularity, he will only see a single edge of the tree and the picture will look like the one 
for a line. With time, he moves away from the tree and more and more intervals corresponding to the edges 
and vertices of the tree will appear. In the limit where his eigentime and his cosmological time go to infinity, 
he will see the image of the whole tree. 

This implies that the observer can reconstruct the domain (up to a global Poincare transformation) from his 
measurements if he waits infinitely long. From his observations, he can reconstruct the edges of the singular 
tree, and — if the spacetime is obtained by Lorentzian grafting on a closed hyperbolic surface — the action of 
the fundamental group on the tree. This amounts to recovering the measured geodesic lamination from which 
the spacetime is constructed. 
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If the spacetime corresponds to a grafted genus g surface and the observer knows the associated Fuchsian 
group (i. e. the hnear part of the holonomy), he can construct the complete domain in finite eigentime. We 
will consider below to what extend the observer can reconstruct the geometry and topology of the spacetime 
without knowing the linear part of the holonomy. 

We now concentrate on the case, where M is a maximal flat globally hyperbolic space-time with closed 
Cauchy surface S of genus g > 2. The universal covering of M is then isometric to a regular domain D. More 
precisely, there is a subgroup G of /som(R^'^) such that D is invariant under the action of G and M = D/G. 
Let r be the subgroup of 5*0+ (2,1) consisting of the SO~^ (2,1) components of elements of G. It is known 
(see [in]) that F is a discrete subgroup of SO'^{2, 1) and that H^/F is a surface diffeomorphic to S. Moreover 
the measured geodesic lamination A dual to the initial singularity of D is invariant under the action of F and 
induces a measured geodesic lamination A on H^/F. 

The main result we present in this section fProposition 15.131) states that if M is not a conformally static 
space-time (which would correspond to the empty lamination) and if its initial singularity is on a simplicial 
tree, then an observer can construct in finite time a finite set of elements of 50"*" (2, 1) which generates a finite 
extension of F. In other words, we will prove the result only when the lamination A is rational (that is, its 
support is a disjoint union of closed curves). We believe that the result could hold also for a general lamination. 
However, in that case some technical issues arise which make the analysis more complex, and we prefer to focus 
on the simpler case where A is rational. 

The main idea is to consider the isotropy group Fq of A 

Fo = {7e^O+(2,l)|7(A)=A}. 

It is clear that Fq is a discrete subgroup of 5'0"'" (2,1) containing F. The quotient H^/Fq is a surface with 
possibly singular points (corresponding to the fact that some element of Fq could fix some point in H^). There 
is a natural projection map 

TT : tf/r ^ eVro 

which is a finite covering. In particular, the index of F in Fq is equal to the cardinality of the fibers of tt and is 
therefore finite. 

Remark 5.7. Any element of Fq is the linear part of an affine transformation that preserves the regular domain 
D. So elements of Fq are the linear parts of the elements in the isotropy group Go of D 

Go = {g e /som(R2'i) | g{D) = D} . 

It should be remarked that in principle there are many subgroups G' of Go (of finite index) such that D/G' is a 
MGH spacetime with compact Cauchy surface. Clearly the intensity of an observer in M is equal to the intensity 
function of some observer in such spacetimes. This suggests that an observer cannot precisely determine the 
group G (or F), but only the group Fo. It should be also noted that in the generic case, F = Fo and it does not 
contain proper cocompact subgroups. 

We say that a leaf / of A is seen by an observer (p, v) if the intersection of D with the support plane orthogonal 
to some point on / intersects I~{p). Note that ii x,y € then the intersection of D with the support plane 
orthogonal to x is equal to the intersection of D with the support plane orthogonal to y. 

In the following proposition (and therefore in the final result of this section) we restrict attention to a 
lamination A with a simplicial dual tree, although it appears quite likely that the proposition holds for general 
laminations. 

Proposition 5.8. Suppose that the dual tree of the lamination A is simplicial. Then the intensity function of 
an observer {p, v) allows one to reconstruct the sublamination A(p.t,) consisting of the leaves of A seen by (p, v) 

Proof. As mentioned above, the intensity function seen by the observer can be split into intensity functions of 
different regions which correspond, respectively, to the edges and to the vertices of the singular tree. In the 
regions corresponding to the edges, the intensity is analytic and behaves as in Example 2 in Section 13.21 It is 
shown there that knowing the intensity on an open subset of is sufficient to determine the positions of the 
edges, and therefore the leaves of the lamination A(p „). □ 

Let us fix a point xq G H^, and denote by Bd the ball in centered at a;o with radius d. For simplicity, let 
us suppose that the point a;o does not lie in a leaf of A. We denote by A^ the sublamination of A made of leaves 
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that intersects Bd'- 

Ad = U / . 
I leaf of A 

Lemma 5.9. For any d > there is a time T such that for t > T the observer {p + tv, v) sees all the leaves in 
Xd, or, equivalently, Xd C Xp+tv.v 

Proof. There is a compact subset K of dD such that ii x ^ Bd then the support plane orthogonal to x intersects 
dD in K. Since {p + tv)r\ D is an increasing sequence of open subsets that cover D, there is a constant T such 
that {p + tv) contains K for t > T. By definition, we then have Ap+t^ C A^, and the conclusion follows. □ 

We now consider the elements of 5*0+ (2, 1) that send leaves of Xd either out of Bd or to other leaves of A^: 

Td = {je S0+i2, 1) I -/{Xd) nBdd Xd} . 
It is easy to check that Tq — C\d>o^d- Note that Td is not discrete. Indeed, transformations 7 such that 
7(Ad) n Bp = form an open subset of 5*0+ (2, 1) that is contained in Td- On the other hand, we will prove 
that the intersection of a neighborhood of the identity with Td is discrete and that this neighborhood can be 
chosen arbitrarily large, by choosing d sufficiently large. 

Lemma 5.10. For any compact neighborhood H of the identity in 5*0^(2, 1), there is a constant d such that 
TdH H is finite. 

Sublemma 5.11. For any a > and d > there is a finite number of strata F of \ X such that F n Bd 
contains a point at distance exactly a from OF. 

Proof. By contradiction, suppose there are countable many strata Fn as in the Lemma, and denote by x„ €E 
Fn n Bd the points such that d{xn, dFn) = a. 

Up to passing to a subsequence, we can suppose that Xn x. li x does not lie in the lamination, then x„ 
definitively lies in the stratum F through x, so Fn ~ F, and this contradicts the assumption on F„. 

If X lies on A, then d{Xn,dFn) = c?(x„,A) — > 0, which contradicts the assumption that this distance is a 
constant larger than 0. □ 

Proof of Lemma \5.1U[ Let do be a fixed number such that Bdo intersects two two leaves h and I2 on the boundary 
of the stratum Fq through xq . 

By the compactness of H, there is a constant r > such that dij2 (a;, 7(3;)) < r for any x G Bd^ and 7 G -ff . 
Note that jlk) intersects Bd with d = do + r. 

If 7 G -ff n Frf with d > do + r, then 7 sends U to some leaves ci and C2 of A^. Clearly, 7(2:0) lies in a stratum 
bounded by ci and C2, and the distance between 7(a;o) and ci is the same as the distance between a;o and li 
(say a > 0). On the other hand, by Sublemma 15. ill there are finitely many strata -Fi, . . . ,F]\f of A^ such that 
Fi n Bd contains a point at distance a from dFi. Moreover, the boundary of each Fi intersects Bd into a finite 
number of segments. 

In particular, there are a finite number of leaves ti , . . . , <m such that every 7 in _ff n sends U to one of the 
leaves t^. However, for two pairs of geodesies (^1,^2) and (ti,t2) in H^, there is at most one isometry sending k 
to ti. This implies that H D Td contains at most 2^'^ elements. □ 

Let us now fix an observer {p,v). It then follows from Lemma 15.91 and Lemma 15.101 that for any compact 
subset H C SO'^{2, 1) and d sufficiently large, there is a time T = T{H, d) such that the observer at proper time 
t > T can list the elements oiTdDH. On the other hand, in principle, an observer would have to wait an infinite 
amount of time to determine if a given element in Tdg lies also in Tq. Indeed, this amounts to determining 
whether such an element lies also in all Td for d > do. The following lemma ensures that this is not the case 
and that the observer can be sure after a finite amount of time that elements oiTdH H also lie in Fq . 

Lemma 5.12. For any compact subset H C 5*0^(2, 1) there is a constant d such that 7 G i? H implies 

Proof. By contradiction, suppose that there is a diverging sequence c?„ and a sequence jn F[ such that 
In e Td„ n H, but 7„ ^ Fq. 

Up to passing to a subsequence, we may suppose that 7^ converges to 700- On the other hand, by Lemma 
15.101 we can choose no big enough so that F^^^ n is a finite set. Now 7„ G F^^ n i? for n > no, and since 
it is a convergent sequence we have that 7„ — 700 for n > ni. This implies that 700 G Fq and contradicts the 
assumption on the sequence. □ 
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We can finally state and prove the main results of this section. 

Proposition 5.13. Let {p,v) be an observer in a domain of dependence D which is the universal cover of a 
Minkowski spacetime obtained by Lorentzian grafting of a closed hyperbolic surface along a rational measured 
lamination. Then the observer can construct a finite set of generators of Fq in finite time. 

Proof. As Fq is finitely generated, there is a compact subset H C 5*0+ (2,1) such that H Fq is a set of 
generators of Fq. By Proposition 15. 81 Lemma fS.lOl and Lemma [5. 121 the observer can detect elements of _ff n Fq 
in a finite time. □ 

5.3. Higher dimensions. In sect ions |5 . 1 1 and we focussed on flat spacetimes of dimension 2 + 1. It appears 
quite possible that an analogous analysis can be done in dimension 3 + 1, or actually in higher dimensions. 
However, proving the results is more involved, since the structure of MGHFC spacetimes is less well understood 
and their description is more complicated than in dimension 2 + 1. 

In dimension 3 + 1, MGHFC spacetimes can still be constructed from a hyperbolic metric on a 3- manifolds 
along with a "geodesic foliation" . Those foliations, however, are different from those occurring on hyperbolic 
surfaces, since they have two-dimensional leaves which can possibly meet along one-dimensional strata. 

Still, it appears plausible that the same conclusions can be reached as in dimension 2 + 1 for observers in a 
domain of dependence which is the universal cover of a MGHFC spacetime in dimension 3 + 1: 

• If the linear part of the holonomy is known, the observer can reconstruct the complete holonomy in 
finite time. 

• The observer can determine in finite time the part of the initial singularity corresponding to the part 
of space he "sees" , which is increasing with time. 

To understand whether those statements hold, it would be necessary to extend to dimension 3 + 1 parts of the 
study done in the previous sections for dimension 2 + 1. 

6. Examples in 2+1 dimensions 

In this section, we study explicitly the light signal reaching an observer in a few simple examples of 2+1- 
dimensional domains of dependence. 

6.1. Explicit holonomies. We consider below three explicit examples, one based on a reflection in the edges 
of a hyperbolic quadrilateral and two based on a hyperbolic punctured torus. In the first punctured torus 
example, the translation part of the holonomy corresponds to a measured lamination with support on a closed 
curve. In the second punctured torus example the support of the measured lamination is more complicated. 

6.1.1. Example 1: A hyperbolic reflection group. In the first example, we consider a group T^/^ which is gen- 
erated by the reflections on the edges of a quadrilateral with angles 7r/3. (This angle condition completely 
determines a presentation of the group, see e.g. [TO]-) The holonomy representation pt : F^/3 — > Isom(M^'^) 
depends on a real parameter t. 

We describe flrst the linear part p\ : F^/3 0(1,2). The construction is based on a quadrilateral p 
with vertices ui,--- ,U4. Consider the hyperbolic plane as a quadric in the Minkowski space R^'^, and let 
wi, W2, twa, ?«4 be the unit spacelike vectors which are orthogonal to the oriented plane through containing 
the geodesic segments (t;i,W2), ("^2,^3), {v^^v^) and (w4,ui). The cosine of the exterior angle of p at Vi is then 
equal to the scalar product between Wi-i and Wi, so that p has interior angles equal to 7r/3 if and only if 
{w^,w^+l) = -1/2 for ah i e Z/4Z. 

The fact that those scalar products are equal means that {wi, W2, W3, W4) is a rhombus in the de Sitter plane. 
In particular it is invariant under the symmetry with respect to a timelike line in M^'^, corresponding to a point 
o G which is the midpoint of both, (wi, v^) and (u2, ^4). 

Choosing a coordinate system compatible with this symmetry, we can write the Wi as 

wi = (sinh(t), cosh(t), 0), ■u;2 — (sinh(i'), 0, cosh(f')), W3 — (sinh(t), — cosh(i), 0), 1^4 = (sinh(t'), 0, — cosh(t')) 

with t,t' satisfying the condition sinh(i) sinh(t') = 1/2. 

One obtains in this manner a representation pt : T 0(2, 1) which can be described as follows. F is generated 
by the elements ai, • • • ,04 corresponding to the reflections in the edges of p, with the relations 

af = 1, (flifli+i)^ = 1 

for all i € Z/4Z. The representation p sends a,; to the reflection in (wi,?;i+i), that is 

p(aj)(x) ^ X - 2{x, Wi)Wi . 
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The quotient of by p(r) is an orbifold. The subgroup r2 of F of elements 7 £ F such that p{'y) is orientation- 
preserving has index two, and the quotient of by p{^2) is a surface. 

There is a unique choice of a deformation cocycle associated to pt, which is obtained by varying t. It can be 
written as Tt = Pf^dpt/dt. So we obtain a one-parameter family of domains of dependence parametrized by t. 
In the following, we mainly consider the simplest case, where t = = sinh~"'^(l/-\/2), so that t' = to. 

6.1.2. Example 2: A punctured torus with a rational measured lamination. Another simple example can be 
constructed, when one chooses as the linear part of the holonomy the holonomy representation of a hyperbolic 
punctured torus. 

The group F is now the free group generated by two elements a, b. We consider the situation with an extra 
symmetry, corresponding to the condition that the images of a, b by the linear part p of the holonomy are 
hyperbolic translations with orthogonal axes. The translation lengths of p{a) and p{b) can then be written as 
2ta,2ti,, subject to the conditions that sinh(ta) sinh(tb) = 1. This corresponds to the condition that the image 
of the commutator of a, b is parabolic. 

In the computations below, we choose, somewhat arbitrarily, the parameters ta = sinh^^(2), tb — sinh^^(l/2). 
We also choose the translation component of the holonomy as the cocycle r corresponding, through the relation 
explained in Section [231 to a closed curve corresponding to 6, with weight 1. 

The domain of dependence obtained in this way is represented on the left in Figure [TJ To compute this image 
— as well as intensity and "distances" to the boundary in other domains of dependence in 2 -I- 1 dimensions, 
figures [5] and [7] below — we use the description of a domain of dependence as an intersection of half-spaces 
bounded by lightlike planes seen in Section [^771 The image is computed by taking a ball of radius 6 in F, for 
the distance defined by the choice of generators described above, and computing for each element of F in this 
ball the corresponding axis and lightlike hyperplanes, and then determining their intersection. 

6.1.3. Example 3: A punctured torus with an irrational lamination. Finally we used a third example in compu- 
tations, which is similar to the previous one in that the linear part of the holonomy is the same. However, the 
translation part of the holonomy is given by a cocycle corresponding to a measured lamination which does not 
have support on a closed curve. The corresponding domain of dependence is given on the right in Figure [TJ 

6.2. Results. Figure [5] presents (in green) the intensity measured by all directions by an observer in the 2 + 1- 
dimensional fiat spacetime described in Section 16.1.11 More precisely, the observer is located above the origin 
at time distances 10, 30 and 50 to achieve the most "readable" results. The Euclidean length of the lightlike 
segments from the observer to the boundary of the domain are drawn in blue. The intensity function becomes 
more complex as the cosmological time of the observer increases. 

Figure [6] is similar, for the 2 -I- 1-dimensional fiat spacetime described in Section 16.1.21 More precisely, the 
observer is located above the origin at time distances 1, 5, 10 and 30. Figure [7| shows the analogous results, for 
the domain of dependence obtained from an irrational lamination, described in Section [6.1.3l 

The computation of both the Euclidean distance and the intensity function are made for an approximation 
of the domain of dependence, as explained above. It follows from Section 4 that the intensity computed in this 
way is not reliable as a continuous function, but only — possibly at least — as a function, due to Theorem 

WM 

Those graphs should be considered as preliminary results, since, even for this relatively simple setting, 
the computations needed to obtain the results are quite involved relative to our programming capabilities 
and computing equipment. It is possible that heavier computations — in particular, computing a better 
approximation of the domain of dependence by using a larger subset of the fundamental group — could lead to 
notably different results. However, it is already apparent in those pictures, that the intensity function behaves 
in a very non-smooth way, as explained in Section 3 and Section 4. 

7. Examples in 3+1 dimensions 

In this section, we consider an example of a domain of dependence in 3+1 dimensions and show that the light 
emitted from the initial singularity and received by an observer contains rich information on its geometry and 
topology. In the first part, we describe the domain of dependence in 3+1 dimensions, based on a construction 
of Apanasov [2] . The second part contains some images of the light emitted by the initial singularity, as seen 
by an observer. 

7.1. An explicit example. We consider a particularly interesting example of domain of dependence, which is 
used in computations below. The remarkable property of this example is that, for only one linear part of the 
holonomy, there is a four-dimensional space of possible translation components. 
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Distances and angles, quadri_pi3, height=10, uu=0, 1=8, npts=200 Distances and angles, quadri_pi3, height=30, uu=0, 1=8, npts=200 

3.5 I ^ ^ 1 1 ^ ^ 1 IS I 1 1 ^ ^ 1 1 1 




Distances and angles, quadri_pi3, height=50, uu=0, 1=8, npts=2Q0 




Figure 5. The intensity (green) and Euclidean distance to the boundary (blue), in a domain 
of dependence based on a quadrilateral with angles 7r/3, with a rational lamination, seen from 
increasing distance from the initial singularity 

7.1.1. The construction of the group. The following example is essentially due to Apanasov 2 . It is a discrete 
group r of Isom{M'^) generated by 8 reflections, so that the quotient H'^/r is a (non-orientable) orbifold of 
finite volume. 

On S"^ = C U {oo} we consider the following circles 

(1) Ck with center at Zk = VSe"^^ and radius 1. 

(2) C with center at and radius 1. 

(3) C" with center at and radius 2. 

It can be shown easily that the configuration of such circles is the one shown in the picture. Moreover the 
angle formed by any two circles in the list (that meet each other) is tt/S. 

We consider the planes Pk,P, P' in H'^ that bound at infinity the circles Ck, C and C". We denote by 7^ the 
reflection along Pk, by 7 the reflection along P and by 7' the reflection along P'. 

Proposition 7.1. The group T generated by 7^,7,7' is a discrete group in Isom{M^) and the quotient H'^/F 
is a non-orientable orbifold of finite volume. 

A fundamental region for the action of F can be obtained as follows. Denote for any k by Ek the exterior of 
the plane Pk, e. g. the region of H'^ \ P that contains 00. Analogously, denote by E the exterior of P, whereas 
/' is the interior region bounded by P' . Then a fundamental region for F is given by 

5 

K ^ f] Ek n E n r . 

k=Q 

As we are interested in R'^'^-valued cocycles, we need to determine explicitly the matrices in 0(3, 1)+ corre- 
sponding to the transformations 7^,7,7'. 
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Distances and angles, torusl, height=l, uu=0, 1=6, npts=200 



Distances and angles, torusl, heiglnt=5, uu=0, 1=6, npts=200 





Distances and angles, torusl, height=10, uu=0, 1=6, npts=200 




Distances and angles, torusl, height=30, uu=0, 1=6, npts=200 




Figure 6. The intensity (green) and Euclidean distance to the boundary (blue), in a domain 
of dependence based on a punctured torus, with a rational lamination, seen from increasing 
distance from the initial singularity 



In R'^'^ we consider coordinates xq,xi,X2,X3 so that the Minkowski metric takes the form —dxg + dx\ + 
dx\ + dx\. Given two real numbers uq, wi and a complex number z = x ~\- iy , vfe denote by (wq, wi, z) the point 
(uo,wi,a;,y) in R^-^ 

We have to fix explicitly an isometry between the half-space model of (denoted by 11 here) and the 
hyperboloid model denoted by H''. Such an isometry (/) : 11 — > H'^ is given by 



(i) = (1,0,0,0) 



d d d 
'>*,iia^ + b— + c-z) = (0,c, a, 6) . 



dx dy d 

With this choice, the plane Pk is given by PkM? C\ v-^, where 

«fe-(3/2,l/2,V3e''^) 
is a unit vector. Analogously, we obtain P = M'^ O and P' = n (w')^ where 

i; = (0,1,0,0) u' = (3/4,5/4,0,0) . 
The associated transformations in 0^(3, 1) then take the form 

jkix) ^ x " 2{x,Vk)vk, j{x) ^ X - 2{x,v)v -/'{x) = x - 2{x,v')v' . 

7.1.2. The construction of the cocycle. In S*^ = CU {oo} we consider the three lines through passing though 
the centers of the circles Ck and denote by W be the union of the planes in H'^ bounding these lines. Clearly, 
W is the union of 6 half-planes which meet along the geodesic Iq joining to oo. We denote these half-planes 
by Wo,Wi,W2,W3,W4,W5 where the indices correspond to the ones of the circles in the obvious way. Let 
Wk e R'^'^ be the vector orthogonal to Wk and pointing towards Wk+i- A direct computation then shows that 
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Distances and angles, torus4, height=l, uu=0, 1=6, npts=200 Distances and angles, torus4, height=5, uu=0, 1=6, npts=200 




Distances and angles, torus4, height=10, uu=0, 1=6, npts=200 Distances and angles, torus4. height=30, uu=0, 1=6, npts=200 




Figure 7. The intensity (green) and Euclidean distance to the boundary (blue), in a domain 
of dependence based on a punctured torus, with an irrational lamination, seen from increasing 
distance from the initial singularity. 



Figure 8. Construction of the linear part of the holonomy 



Wk is given by 

Wk = (0,0,^6^^^). 

Note that Wk+s — —Wk, where the index k is considered mod6. 

Now the F-orbit of VF is a branched-surface in H''. In particular, the sets W = T ■ W and Wq = F • have 
the following properties: 
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• Wo is a disjoint union of geodesies. 

• Every connected component ofW\ Wq is a convex polygon (with infinitely many edges) and every edge 
is an element of Wq. 

• There are exactly six faces up to the action of F. Indeed, let Fk be the face of W bounding and 
contained in Wk- Then the orbits oi Fq, . . . , are disjoint and cover W. 

• Every connected component of H'^ \ is a convex polyhedron. 

These properties can be proved by considering the intersection of W with K and using the fact that W is 
orthogonal to the faces of K. 

We can then use a general construction explained in [7^ to obtain non-trivial cocycles. Given six numbers 
such that 

5 

(39) = 0, 

i=Q 

we obtain a cocycle via the following prescription. We associate to the face Fi the number a^. In this way a 
number a{F) is associated to every face F by requiring that a{a{F)) — a{F) for every a £ T. 

Then we fix a basepoint xq in H'^ that does not lie in W. Given a transformation a e F we construct a vector 
in M"^'^ in the following way: we take any path c joining xq to a{xo) and avoiding Wq. The path c intersects 
some faces F^, . . . We consider the unit vector S M'^'-^ orthogonal to F^ and pointing towards a{xo) 

and set 

n 

r(a) = ^a(i^^')w^' . 
j=i 

It can be easily checked that 

• T(a) does not depend on the path c (this essentially follows from (P^]) '). 

• T is a M'^' ^-valued cocycle 

• changing the basepoint changes r by a coboundary. 

Let H C M.^ he the subspace of solutions of ([55]) . which is of dimension diuiH = 4. From [2, we have that the 
map 

H i?^(r,M^'^) 

is injective. This map can be computed explicitly. More precisely, we fix the base point xq in the region of K 
between Wq and Wi. 

Given a set of numbers oq, oi, 02, 03, 04, 05 that satisfy pop we compute the corresponding cocycle r evaluated 
on the generators. This yields 

r(7) = r(7') = t{jo) = ^-(71) = 

r(72) = wi - 72W1 = 3ait;2 

r(73) = wi - J3W1 +W2 - 73102 = 3(ai -I- 02)^3 

'^(74) = -Wo + l4,wo -W5+ 74W5 = -3(ao + a5)w4 

'^(75) = + 75 Wo = -3aoW5 ■ 
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7.2. Computations. The computations of ttie intensity functions were limited by the speed of the available 
computers, and more complete computations could provide better results. As in dimension 2 + 1, we constructed 
a domain of dependence in M^'^ , invariant under the group actions described above via the construction in Section 
12.71 However, the computations are much heavier in dimension 3 + 1, so we only considered the elements of the 
fundamental group in a ball of radius 4. 

Although we only considered one linear part of the holonomy — the one in Section 17.1.11 — we worked with 
two deformation cocycles, one corresponding to weights (1,0, 0, 1, 0, 0) as described in Section r7.1.21 the other to 
the weights (1, 1/2, 0, 1, 1/2, 0). In both cases, the observer was located at the point of coordinates (50, 0, 0, 0). 
This is a somewhat arbitrary choice, made after trying different possibilities, which leads to interesting pictures. 

The intensity measured by the observer for the first choice of cocycle is presented in Figure [TUl with different 
colors encoding different values of the intensity. 
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Figure 10. Computed intensity, with translation coefficients (1, 0, 0, 1, 0, 0), from an observer 
at an increasing distance from the initial singularity 

It should be noted that those results are even less certain that those obtained in dimension 2 + 1 above. This 
is due to the fact that we compute the limit intensity for a decreasing sequence of finite domains of dependence 
approximating the domain under examination. In dimension 2 + 1, Theorem 14.131 and Proposition l4.22l indicate 
that the limit intensity is the intensity of the limit, if the limit is the universal cover of a MGHFC spacetime. 
By contrast, in dimension 3 + 1, we only know by Theorem 14. 141 that the intensity of the limit is at least equal 
to the limit intensity, and at most equal to three time the limit intensity. So the intensity functions computed 
here are the limit intensity (which is a well-defined notion for any domain of dependence, see Theorem I4.14p 
which differs from the "real" intensity by a factor at most three. 

The limit intensity computed for the second choice of cocycle is in Figure 1111 It is apparent how the less 
symmetric cocycle leads to a distortion in the figure. The symmetry of degree six, which is present in the linear 
part of the holonomy, is readily apparent in Figure [TUl In Figure [TT] it remains visible, but with differences in 
the size of the corresponding parts of the picture. 
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Figure 11. Computed intensity, with translation coefficients (1, 1/2, 0, 1, 1/2, 0), from an ob- 
server at an increasing distance from the initial singularity 



Even for this fairly simple example, it would be interesting to perform more powerful and complete compu- 
tations, for instance by computing the domain of dependence with all elements of the fundamental group in a 
ball of radius larger than 4. It is conceivable that one would obtain somewhat different pictures. Additionally, 
the picture should vary with the position of the observer. It should be simpler for an observer close to the initial 
singularity, but become increasingly complex as the observer moves away from it. 
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